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ABSTRACT 

The  purpose  of  this  investigation  was  to  study,  by  way  of 
permeation  experiments,  structural  differences  between  synthetic  poly- 
els  1,4-isoprene  (SNR)  and  natural  rubber  (NR),  as  well  as  differences 
resulting  from  filling  with  carbon  black.  The  study  aimed  further  at 
testing  Frisch's  theory  of  isotope  effects  in  activated  diffusion  and  on 
the  unit  size  of  the  free  volume  in  rubbers. 

The  rubber  samples  consisted  of  thin  croasllnked  sheets  of  NR 
and  SNR.  The  carbon  content  (N330  HAF  black  of  about  300X  particles 
diameter)  in  the  filled  samples  was  16%  by  volume. 

The  permeation  rate  of  and  Dj  through  filled  and  unfilled 
NR  and  SNR  were  measured  by  the  time-lag  method  as  a  function  of  tempera¬ 
ture  and  pressure  (AO-Sj^C  and  1-55  cm  of  Hg).  The  permeability  coefficient 
P,  diffusion  coefficient,  D,  and  solubility,  S,  thus  deduced,  were  further 
used  to  evaluate  the  energies  of  activation  for  permeation  and  for  diffusion 

E  and  E.,  and  the  heat  of  solution,  AH  according  *o  the  Arrhenius 
pa  8 

equation. 

The  results  show  that  has  generally  a  higher  permeability 
and  diffuslvity  than  D^,  and  that  the  permeabilities  or  diffusivlties 
are  markedly  reduced  by  the  presence  of  the  filler.  On  the  other  hand, 
the  solubilities  in  the  filled  rubbers  were-  found  to  be  almost  double 


that  of  the  unfilled  samples. 

The  diffusivities  and  the  activation  energies  for  diffusion 
and  permeation  are  largely  Independent  of  the  gas  pressure.  The  heats 
of  solution,  *,H#,  for  both  gases  in  the  unfilled  samples  are  all  positive, 
but  in  the  filled  samples,  AH#  is  positive  for  D2*  and  negative  for  Hg. 
These  signs  of  A.Hg  are  consistent  with  their  solubilities. 

A  qualitative  analysis  of  the  heat  of  solution  permits  one  to 
estimate  the  heat  of  adsorption  on  NR  to  be  approximately  800  cal/mole; 
the  overall  heat  of  solution  further  allows  the  estimation  of  2  lccal/mole 
for  the  energy  of  hole  formation  in  the  rubber. 

The  solubilities  (adsorption)  of  Hj  and  Dj  were  found  to  be 
low  in  both  rubbers  and  on  the  carbon.  The  adsorption  of  on  the  rubber 
is  about  3.2  x  10  6  g./g.  of  rubber,  and  8  x  10~*  g./g.  of  carbon.  The 
latter  is  probably  completely  wetted  by  rubber. 

Based  on  the  quantum  mechanical  treatment  of  diffusion  suggested 
by  H.  L.  Frisch,  the  unit  free  volume  in  the  unfilled  rubbers  is  calcu- 
lated  from  the  value  of  6AHg,  i.e.,  the  difference  of  AKg  for  D2  and 
in  the  same  medium  to  be  51  and  64  X3  respectively  for  NR  and  SNR. 

However,  Frisch's  equation  can  not  be  applied  to  calculate  the  free  volumes 
in  the  filled  staples  since  AHg  was  found  to  be  positive  for  this  case. 

The  jump  length,  X,  of  l^in  NR  and  in  SNR  are  found,  according 
to  the  absolute  rate  theory,  to  be  almost  the  same  (about  21  X)  in  both 
rubbers,  whereas,  X,  for  about  26  A°  in  NR  and  30  X  in  SNR. 

A  semi-quantitative  analysis  of  the  diffusion  mechanism  in  the 
filled  samples  shows  that  the  data  of  P,  D,  and  S  for  the  filled  samples 


li 


can  be  related  to  those  of  the  unfilled  samples  accounting  for  the 
reduction  of  the  cross  section  of  the  filler,  the  imposed  tortuosity 
of  the  diffusion  path,  and  the  increase  in  the  viscosity  of  the  rubber 
due  to  the  presence  of  the  filler  particles. 
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I.  INTRODUCTION 

A.  Introduction  and  Purpose  of  the  Thesis 

The  permeation  of  small  molecules  through  a  polymer 
is  a  process  of  mass  transport.  The  measurement  of  permeabi¬ 
lity  may  be  utilized  to  characterize  the  molecular  or  .physical 
structure  of  a  polymer  and  even  reflect  on  mechanical  responses e 

The  mechanism  of  the  movement  of  a  gas  through  a 
polymer  is  usually  conceived  as  an  initial  sorption  and  disso¬ 
lution  of  the.  gas  in  a  very  thin  surface  layer,  followed  by 
diffusion  of  gas, molecules  under  a  concentration  gradient,  and 
eventually  evaporation  into  the  gaseous  space  at  the  other  sur¬ 
face.  In  the  case  of  vapors  the  step  prior  to  the  dissolution 
will  be  a  process  of  condensation  rather  than  of  adsorption. 

The  situation  will  then  be  complicated  by  such  facts  asi  their 
solubilities  are  not  directly  proportional  to  the  pressures  and 
the  diffusion  coefficients  often  become  highly  dependent, on  the 
concentration  of  the  penetrants  in  the  polymer.  Furthermore,  vapors 
in  the  polymer  also  act  as  plasticizers  or  swelling  agents, 
lower  the  glass  transition  temperature  and  loosen  the  segmental 
.motions  of  the  polymer. 

By  contrast,  the  situation  is  usually  much  less  com¬ 
plicated  for  the  permeation  of  non-condensing  gases,  in  most 
cases  Henry's  law  is  obeyed  and  the  diffusion  coefficients  are 
independent  of  the  concentration.  For  these  reasons  gases  are 
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preferred  for  studies  into  the  mechanism  of  transport  via 
sorption  and  diffusion,  and  for  the  information  derived  there¬ 
from  on  polymer  structure,  morphology,  and  othtt*  features*  Although 
information  obtained  from  the  studies  of  the  effect  of  polymer- 
penetrant  interaction  on  the  transport  process,  using  vapors, 
is  of  great  value  to  general  membrane  studies,  i.t  will  be 
treated  here  only  tangentialiy. 

The  quantities  obtainable  from  gas  transport  studies 
are,  depending  on  the  theories  of  mechanism  applied,  permea¬ 
bility,  diffusivity,  solubility,  activation  energies,  jump 
» 

distance  and  the  free  volume  in  the  polymer.  The  magnitude  of 
these  quantities  characterizes  the  penetrant-polymer  system 
and  is  affected  by  the  polarity  or  size  of  the  penetrant  and 
the  nature  of  the  polymer  (i.e.,  polar  or  side  groups  along 
the  chain,  double  bonds  in  the  chain  backbone,  degree  of  crcss- 
linking,  degree  of  crystallinity,  plasticizers,  and  fillers). 
Any  change  in  the  morphology  or  structure  of  the  polymer  will 
cause  changes  in  these  quantities. 

It  is  the  purpose  of  this  study  to  evaluate  all 
these  quantities  for  four  rubber  samples,  namely  filled- 
and  unfilled  natural  rubber  and  synthetic  poly-cis- 
isoprene,  as  a  function  of  pressure  and  temperature ,  using  K2 
gas  ar.d  its  isotope,  D2,»  as  penetrants.  The  isotope  effect 
that  will  be  shown  by  the  permeation  data  will  then  be  utilized 
to  evaluate  the  free  volumes^1),  and  the  effects  of  p  and  T  on 
the  penetration  rates  will  be  evaluated  in  terms  of  the  other 
quantities  mentioned.  All  rubber  samples  were  prepareu  under 
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the  same  conditions  and  have  practically  the  same  chemical 
composition,  i.e.,  poly(cis  1 ,4-isoprene) ,  but  were  selected 
because  synthetic  and  natural  rubber  have  been  found  to  show 
significant  differences  in  their  tensile  strength,  elongation 
and  modulus.  The  filled  rubbers  contained  carbon  black  and 
were  designed  for  study  of  thp  effect  of  the  filler  on  the 
transport  process. 

The  information  gained  from  this  study  will  be  in¬ 
terpreted  in  terms  of  rubber  structure  and  rubber-filler 
interaction  and  should  be  of  obvious  theoretical  and  practical 
importance, 

B,  Basic  V'cuntions  of  Permeability  of  he*es  through  Polymeric 
Membranes 

According  to  Fick’s  first  law  of  diffusion,  mass 
transport  or  the  flow  rate,  Q  in  ml“2t"1,  of  a  substance  in 
a  medium  is  proportional  to  its  gradient  of  concentration* 

Q  =  -D  ~  (1) 

dx 

The  proportionality  factor,  D,  is  called  the  diffusion  coeffi¬ 
cient*  However,  this  equation  gives  only  the  steady-state 
condition  for  diffusion  in  dilute  systems.  If  the  concentra¬ 
tion  changes  with  time  as,  for  instance,  during  the  transient 
period  of  flow,  the  formalism  of  the  equation  has  to  be  changed. 
Consider  a  region  of  unit  cross-section  and  of  length 
dx,  evtending  from  x  to  x+rtx,  as  shown  in  Fig,  1, 


X 


Fig.  1  Diagram  of  transient  flow 


The  increase  in  concentration  within  this  region  in  unit  tine 
divided  by  the  volume  dx,  is  the  excess  of  material  diffusing 
into  the  region  over  that  diffusing  out, 
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therefore, 


(2) 


This  equation  is  Fick's  second  law  of  diffusion  in  one  dimen¬ 
sion.  For  three  dimensions,  it  takes  the  form  of 


is. .  D(  i2|  .  £§  +  ils.  j 

^t  '  *x2  *  y2  j)  a2 


D72c 


(3) 


providing  D  is  a  constant.  However,  in  many  systems  D  depends 
markedly  on  the  concentration.  In  that  case,  one  has  to  write 
equations  (2)  and  {3)  as 


where  D  is  now  a  function  of  x,  y,  2  and  c.  The  aho^o 
tions  to  he  solved  explicitly,  under  th*  expert *i»onta1 
conditions,  to  obtain  D. 

Co  Methods  for  Solutions  to  the  Kiustions 

1.  Stcadv^tste  Condition 

Consider  a  substance  flowing  through  a  membrane, cf 
thickness, L#  where  the  concentrations  and  c2  of  the  nubc- 
tance  at  the  inlet  and  outlet  surfaces  >;f  the  membrane  are 
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kept  constant.  In  the  steady  state,  the  concentrations  of 

all  sections  of  the  membrane  remain  unchanged,  provided  the 
diffusion  coefficient,  D,  is  constant.  The  diffusion  equation 

(2)  becomes  then 


*c  _  d2c 
3t  D  dx2 


0  , 


or 


(6) 


On  integrating  the  equation  twice  with  respect  to  x,  and 
introducing  the  conditions  that  at  x  *  0,  c  *  clt  and  at  x  *  l» 
c  *  c2,  one  obtains 


c  -  c. 
- i  *  ^ 

e2-  ci  L 


(?) 


Clearly,  the  concentration  changes  linearly  from  c1  to  c2 
through  the  membrane  as  long  as  the  membrane  is  uniform. 

Since  the  flow  rate  of  the  substance  is  constant  at 
the  steady-state ,  equation  (1)  can  also  be  integrated  to  give 


P(ci  -  c2) 


as  long  as  D  remains  constant.  If  the  thickness,  L,  and  the 
surface  concentrations,  ci  and  c2,  are  known,  D  can  be  determined 
from  a  single  observation  of  the  flow  rate. 

If  a  gas  or  vapour  diffuses  through  a  membrane,  and 
the  pressures  Pi  and  p2  at  the  two  sides  of  the  membrane  are 
known,  the  steady-state  flow,  Qs»  is  given  by 


? 


p  (Pi  -  P2) 


(9) 


where  P  is  the  proportionality  constant  termed  permeability, 
with  units  of  cm  sec*1  and  Qg  is  in  units  of  gm  cm“2see“1. 
Since  the  gas  rarely  dissolves  to  any  extent  in  a  membrane  and 
there  is  a  linear  relationship  between  the  external  gas  pressure 
and  the  corresponding  concentrations  within  the  surface  layer 
of  the  membrane,  it  is  possible  to  express  the  surface  concen¬ 
tration  by  Henry's  law, 


c  *  Sp  (in) 

where  S  is  the  solubility  in  cm"1.  Substituting  the  two  sur¬ 
face  concentrations,  c^  *  Sp^,  and  c2  *  Sp2,  into  equation  (9) 

and  comparing  with  equation  (p ) ,  one  readily  obtains  the 
relation, 

P  *  DS  (11) 


which  indicates  that  the  transmission  rate  is  the 
result  of  two  concurrent  processes,  namely,  dissolution  and. 
diffusion,  and  depends  on  the  product  of  the  factors  describing 
these  processes. 

New  let  us  consider  the  case  where  the  diffusion 
coefficient  varies  with  the  concentration.  Fick's  first  law 
may  still  be  considered  to  hold  if  one  permits  D  to  be  a  mean 

value  deduced  from  the  measurement  of  the  steady-state  flow 
rate.  Integrating  the  equation  between  the  two  surface  con¬ 
centration,  ci  and  Cg,  yields  then 
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s 


'i  / 


2  D  (C4-C9) 

Ddc  = - - — — 


(1?) 


where 


D  ■ 


cl~c2  ^ 

Co 


Ddc 


(I'M 


2*  Time  Lag 


The  most  general  solution  to  Pick's  second  law  is 


given  by  a  Fourier  series.  Daynes^2^  and  Barrerw/  showed, 
under  the  correct  boundary  conditions  and  with  a  constant  or 
average  D,  that  the  amount  of  diffusant,  in  units  of  ml"2, 
which  passes  through  a  membrane  in  time,  t,  is  given  by 


.(3) 


*5£i.ft  -  IlL)-  2**pl  y  C”1)11  /-Dn27T2t 

T'  to'  »rnt'1V 


/- Dn  7 T  t  \ 

expi— p— ) 


(it) 


When  t  is  sufficiently  large,  the  exponential  term  becomes 
negligibly  small  and  the  equation  reduces  to 


(15) 


A  plot  of  Qt  versus  t  will  result  in  a  straight  line  at  the 
longer  times  which  represents  the  regime  of  steady  state  flow. 
The  slope  of  this  gives  the  constant  flow,  Qg  in  units  of 
ml^t"1  since  L  is  a  constant.  Extrapolating  this  line  to 
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the  t-axis,  the  intercepts  ,  is  given  by 


from  which  D  can  be  calculated.  Once  the  two  quantities  D 
and  P  are  known*  the  solubility  S  follows  readily  from  equa¬ 
tion  (11).  This  is  a  remarkable  method  to  evaluate  all  three 
essential  quantities. 

As  mentioned*  the  foregoing  derivation  requires  a 
constant  D.  This  is  probably  the  case  for  the  diffusion  of 
small  molecules  in  a  polymer.  For  organic  vapors*  or  water 
vapour*  the  diffusion  constants  are  often  dependent  on  the 
concentration  of  the  penetrant  in  the  membrane.  This  means 
that  the  diffusivities  can  not  be  readily  obtained  from  this 
so  called  time-lag  method. 

When  D  depends  on  concentration,  and  c  *  c(x,t) 
denotes  the  concentration  of  gas  in  the  membrane  at  a  distance 
x  from  the  membrane  surface,  on  the  high  pressure  side,  at 
time  t,  Frisch^)  has  shown,  by  solving  equation  (4)  under  the 
conditions! 

c(x,0)  ■  0  for  x>0 
c(0,t)  ■  c0  for  t»0 
c(Lit)  *  0  for  t>0 

where  cQ  is  the  concentration  of  gas  on  the  high  pressure  side* 
that  the  time  lag,  t,  is  given  by 
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r 


A. 

/ 


x  cs(x)dx 


/ 


D(u)du 


(1?) 


The  concentration  inside  the  membrane  at  the  steady  state, 
cg(x),  can  be  found  in  principle  as  a  solution  of  the 
quadrature 


D(u)du  «  - 
L 


/■ 


D(u)du 


(18) 


where  u  is  a  dummy  variable,  and  the  relation  of  the  diffusion 
coefficient  and  of  the  concentration  must  be  of  a  known  form, 
or  bo  assumed  to  suit  the  analysis# 

The  use  of  equations  (17)  and  (18)  can  be  illustrated 
for  two  cases i 

In  the  first  and  simpler  case,  where  D  is  constant, 
i*e#,  D(c)  ■  D0  *  constant,  equation  (18)  becomes 

cs(x)  *  c0  (1-  |  ) 

When  this  is  substituted  into  equation  (17)  with 


o 

D(u)du 


D  qCq 
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one  obtains  the  known  result 


r  *  L2  /  60o 

For  the  second  case,  where  D(c)  *  D0  (  1+bc),  D0  is 
the  diffusivity  at  zero  concentration  and  b  is  a  constant, 
Frisch  has  evaluated  c8(x)  from  equation  (18)  through  a  series 
of  expansion  in  powers  of  b.  For  small  values  of  b,  an  appro- 
ximate  expression  for  f  was  obtained, 


60o\l  +  »bc0  / 

Obviously,  DQ  and  b  must  be  determined  by  measuring  the  time 
lag  at  two  concentrations, 

3*  Short-Time  Approximation 

When  the  diffusion  coefficient  is  small,  it  would 
require  an  inconveniently  long  time  to  establish  a  steady 
state  of  permeation.  In  this  case  the  diffusivity  may  be 
obtained  from  the  intercept  of  the  versus  t  plot. 

An  alternative  solution  to  the  diffusion  equation 
(2)  at  short  times,  for  constant  D  has  been  proposed  as (5) 


£P 

dt 


2APjS  /  D  \* 

v  lirtl 


(1°) 
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where  is  the  rate  of  increase  of  pressure  on  account  of 
at 

the  diffusant  transport  through  a  membrane  of  area,  A,  into 
an  initially  evacuated  space  of  volume  vj  S  is  the  solubility 
and  p^  is  the  gas  pressure  at  the  inlet  surface  of  the  mem¬ 
brane.  When  t  is  sufficiently  small,  only  the  first  term  of 
the  summation  in  equation  (19)  is  important 


A  plot  of  lnjt^  j  against  l/t  should  be  a  straight  line  of 

slope  f-  £-1  from  which  the  value  of  D  is  obtained.  This 
method  was  used  by  Rogers,  Burity  and  Alpert'-''  and  has  been 


applied  successfully  in  several  quite  different  systems  by 
Meares^  • 


at ion  and  Desorption 


The  diffusion  coefficients  of  vapours  in  polymers 
are  several  orders  of  magnitude  lower  then  those  of  gases. 

The  estimation  of  a  steady  rate  of  permeation  is  therefore 
very  time  consuming.  However,  the  kinetic  study  of  the 
sorption  and  desorption  of  gases  or  vapours  in  polymers  can 
be  made  into  a  simple  and  rapid  method  for  the  determination 
of  diffusion  coefficients. 

The  typical  apparatus  employed  by  Garret  and  Park^) 
is  shown  in  Fig.  2.  A  specimen  of  known  shape  is  suspended 
from  a  quartz  spiral  balance  in  an  atmosphere  of  the  gas  or 
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vapour  to  "be  studied  at  constant  pressure  and  temperature. 
The  rate  of  increase  in  weight  of  the  sample  is  determined 
from  the  extension  of  the  spring  as  a  function  of  time. 


Fig.  2  Quartz  spring  sorption  apparatus. 

A,  Stirred  vapour  source.  B,  magnetic 
stirrer.  C,  cold  thimble  for  desorption 
studies.  D,  quartz  spring,  S,  polymer 
specimen  on  glass  cradle, ~  F,  ball  joint 
to  vacuum  system. 


When  D  is  independent  of  concentration,  the  amount  of  diffu- 

sant,  Q  ,  taken  up  by  a  planar  sheet  in  time,  t,  is^ 
t 

q  g  exp  j-Tuk+l)  IT/  l]  Dt  j 

— iw -  (2,)  ■ 

k*0 


14 


where  Q«  is  the  equilibrium  sorption  attained,  theoretically, 
after  infinite  time.  The  same  equation  holds  also  for  desorp¬ 
tion  if  Qw  is  the  amount  desorbed  after  infinite  time  and 
is  the  amount  lost  up  to  time,  t.  For  long  times,  equation 
(21)  may  be  approximated  by 


For  short  times  the  approximation 


(23) 


is  adequate,  A  plot  of  Qt/Q»  against  (t/L2)^  is  initially 
linear  so  that  D  can  be  evaluated  from  the  slope.  Also  if 
ti  is  the  time  at  which  Q^/Q*  *  then 


L2  c  4,919  x  10"2 

64H  (  ) 


(24) 


When  D  is  a  function  of  concentration,  the  quantita¬ 
tive  relationship  between  the  diffusion  coefficient  and  the 
concentration  can  be  determined  from  a  series  of  sorption 
experiments  at  different  initial  concentrations.  For  each 
sorption,  the  Q^/Q.  versus  (  t/L2)^  plot,  according  to  equation 
(23),  gives  3ome  mean  value,  D,  of  the  variable  diffusion  coeffi 
cient  averaged  over  the  range  of  concentrations,  in  the  membranes 
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Calculations  have  shown  that,  for  any  one  experiment  a  reason¬ 
able  approximation  for  D  is  given  by 

Ddc 


Using  this  relationship  a  graph  of  Dc0  as  a  function  of  c0 
can  be  drawn.  Numerical  or  graphical  differentiation  with 
respect  to  cQ  gives  a  first  approximation  to  the  relationship 
between  D  and  cv  In  many  cases  this  first  approximation 

is  sufficiently  accurate.  But,  it  is  sometimes  possible  to 
obtain  a  better  approximation  by  the  simpler  method  of  using 
both  sorption  and  desorption  data.  If  Da  and  Dd  are  calculated 
from  the  initial  gradients  for  sorption  and  desorption,  respec¬ 
tively,  then  *(DS  +  Dd)  is  a  better  approximation  to 

r° 

(  l/c0  )  I  Ddc 
0 

-  (9) 

than  either  Ds  or  Dd  • 


Temperature  Dependence 


The  rate  of  transport  of  geses  or  vapours  in  a  poly¬ 
mer  varies  with  temperature.  In  most  cases,  the  temperature 
dependence  of  P,D»  and  S  follows  the  Arrhenius  equation  as  a 
reasonably  good  approximation.  They  can  then  be  represented  by 


_  .  -Ep/RT 

P  ■  P  e  * 
o 

D  *  D0  e“E^T 

S  *  S  e"®s/RT  m  g  0-  AHt./r;T 
o  o 


(25) 


where  P0,  DQ  and  SQ  are  constants j  E  .  ?.u  -  '  ~£Q  are  experir 
mentally  determined  activation  energies:  for-  iv.eation,  diffu* 


sion  and  dissolution.  The  activation  em-r-y  of  dissolution 
is  customarily  referred  to  as  heat  of  AHS*  The 

above  equations  may  be  written  in  terms  of  th*  differentials, 
obtained  from  the  slopes  of  the  semi-lp.-^ri-  *..> c  graphs,  as 


-Rd  In  P 


-R  d  l  h  .  D 


(26) 


-R  d  In  S 
d  (1/T) 


The  heat  of  solution  can  be  readily  de  teruinr  >.  from 


AHS  *  Ep  -  Ed 


(27) 


since  P  *  DS. 


Thus,  the  various  energies  can  be  obtained  from  the 
slope  of  the  plots  of  the  correspond;:  n,;  1<>  .•  '  thmic  quantities 
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versus  l/T.  The  activation  energy  for  diffusion  is  associated 
with  different  meaning  in  various  diffusion  theories,  as  shall 
be  discused  in  the  next  section. 

E.  Theories  for  Diffusion 

There  are  several  theories  dealing  with  the  problem 
of  transport  of  gas  through  polymer  membranes  and  the  role  of 
the  activation  energy  of  diffusion.  The  volume  fluctuation 
theory  postulates  that  diffusional  jumps  of  the  penetrants 
arise  from  volume  fluctuations  and  that  the  experimental 
activation  energy  is  an  average1 value  for  all  fluctuations 
large  enough  to  allow  a  jump^10),  The  hole  formation  theory 
assumes  that  the  polymer  molecules  have  to  be  pushed  far 
enough  apart  to  form  a  hole  sufficient  for  the  infiltration 
of  the  gas  molecules,  the  activation  energy  in  this  theory 
is  the  energy  needed  to  move  the  chains  apart.  The  transition 
state  theory  postulates  that  the  diffusants  must  be  thermally 
activated  in  order  to  surmount  a  potential  energy  barrier 
encountered  in  the  course  of  their  jumps  from  one  equilibrium 
position  within  the  membrane  to  the  next.  Since  the  interpre¬ 
tation  of  our  data  is  based  on  the  transition  state  theory, 
we  shall  examine  it  in  some  detail. 

The  transition  state  theory  as  developed  by  Eyr^ng 
and  co-workers^ *3)  states  that,  in  diffusing  in  solution 
(solid  or  liquid  ),  a  molecule  of  solute  and  one  of  solvent 
are  required  to  slip  past  on''  ane+b»r.  the 

distune*  between  two  successive  equilibrium  positions  is  X  . 
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The  change  of  the  standard  free  energy  with  distance  can  then 
be  represented  by  the  curve  in  Pig.  3. 

Directi.on  of  Diffusion 
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Pig.  3  Energy  profile  for  diffusion 


For  ideal  or  highly  dilute  solutions  at  rest,  the  standard 
free  energy  is  the  same  in  all  equilibrium  positions.  Thus  the 
energy  of  activation  will  be  the  same  in  the  forward  and  back¬ 
ward  direction,  the  specific  rate  constant,  w,  is  the  same  for 
transport  in  either  direction,  and  the  molecules  move  at  random 
with  no  net  transport.  If  there  is  a  concentration  gradient, 
dc/dx,  and  X  is  the  distance  between  the  successive  equilibrium 
positions,  the  concentration  in  moles  per  cubic  centimerof  solute 
in  the  initial  and  final  states  of  diffusion  jumps,  are  c  and 
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c+  ~X  respectively,  The  number  of  molea  of  solute  moving 
in  forward  direction,  e,g.,  from  left  to  right,  through  a  unit 
cross  sectional  area  per  second  is  then  given  by 

-2  -1 

Qf  B  Vf  *Xow  moles-cm  *-sec  A 

where  w,  expressed  in  terms  of  concentrations,  is  the  specific 
reaction  rate  for  diffusion,  i.e.,  the  number  of  moles  of  solute 
molecules  moving  from  one  position  to  the  next  per  unit  time. 
Similar.ily,  the  rate  of  movement  in  the  backward  direction,  is 

Qb  3  vb  “  X(c+X~)w  mole -cm“2 -sec*"1 

The  resultant  flow  from  left  to  right  i3  therefore 

Q=V**«-X2w~  moles-cm*’2-sec*’1 

On  applying  Fick's  first  law, 

Q  *  -  D  ~  moles-cm*,2-sec"1 


where  k  and  h  are  the  Boltzmann  and  Planck  constants,  and  K* 
is  the  concentration  equilibrium  constant  for  the  formation 
of  the  activated  species* 

Introducing  the  standard  free  energy 
AF°*»  -RT  In 

it  follows  that 

w  „  kT  e-  AF°*/RT  b 
h 

since  AF0*  *  AH0*  -  T  AS^ 

« 

chosen  to  be  the  state  of  unit  concentration  because  rate 
constants  are  usually  expressed  in  terms  of  concentrations* 

The  quantities  AF°  ,  AS0*  and  AH0*  are  called  the  free 
energy,  the  entropy  and  the  enthalpy  of  activation  respectively. 

The  experimental  activation  energy  for  diffusion, 

Ed,  is  obtained  from  the  Arrhenius  equation, 

d  In  w 

dT  RT2 


K, 


♦  J 

kT  AS^/R  -AH°/RT 

t  e  e 


(29) 


•  The  standard  state  may  be 


Comparison  with  equation  (29)  shows  that 
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Accordingly,  from  equations  (28),  (29)#  and  (30)#  the  diffu¬ 
sion  coefficient  may  be  expressed  by 

D  -  .  X2  U  0AS°4/R  e-E^T  Ol) 

h 

The  rate  constant,  w,  can  be  written  in  terms  of  the  partition 
functions  of  the  system  in  the  activated  and  normal  states  at 
0°K,  and  hence  equation  (28) then  becomes 

D  -  X  2w  -  X2  2*e‘  *oAT  {32) 

h  P 

from  which  the  absolute  magnitudes  of  diffusion  coefficients 
can  in  principle  be  determined  by  evaluation  of  the  partition 
function^  ^  • 

If  the  degrees  of  the  freedom  of  the  diffusional 
movement  are  assumed  to  be  translational  only,  i»e«,  if 
the  rotational  and  vibrational  contributions  in  the  activated 
and  normal  states  are  taken  to  equal,  the  translational 
partition  function  for  three  degrees  of  freedom  within  the 
free  volume,  Vf,  is  given  by 

F  ./  2irmkT  '3/2 

tr  I  ^3^  )  f 

The  activated  species  has  only  two  degrees  of  freedom  so  that 
its  translational  partition  function  is 
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Setting  Ed=N?0  and  substituting  for  the  ratio  fVv in 
equation  (32)#  one  finds  for  D 


(33) 


Comparing  equations  (31)  and  (33)  leads  to 


v  1/3 
f 


h  e 


-(  £S0+/R+l)/(27rmkT)* 


(34) 


This  equation  can  be  used  to  calculate  the  standard  free 
entropy  of  activation,  if  the  free  volume,  Vf,  is  known. 
Frisch  has  developed  an  innovative  method  for  the  estimation 
of  the  free  volume  base  on  the  quantum  isotope  effect  in  per¬ 
meation.  This  is  discussed  in  the  following  section. 

F.  Quantum  Isotope  Effect  in  Permeation 

Hydrogen  and  its  isotope,  deuterium  have  the  same 
effective  molecular  volumes  notwithstanding  the  difference 

in  their  mass.  When  such  small  gas  molecules  are  dissolved 
in  microscopic  cavities  in  the  polymer  matrix,  they  may  be 
thought  of  as  free  particles  confined  in  a  potential  energy 
box  in  the  quantum  mechanical  sense.  If  the  .cavity  is  approx¬ 
imated  by  a  cube,  the  quantized  translational  energy  for 
the  particle  is  given  (in  three  dimensions)  by 

3  n2h2 

* n>i  * 

3  h2 

,  *  - - 5  <  35) 

i  8  mjd^ 


or 


€ 
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At  temperatures  not  far  from  room  temperature  the  particle  is 
in  its  ground  state with  the  principal  quantum  number, 
n  «  1.  Equation  (35)  thus  relates  the  dimension  of  the  cavity, 
d,  and  the  mass,  m^,  to  the  potential  energy,  c^»  of  the  parti¬ 
cle,  Realizing  this,  Frisch^  has  proposed  a  model  based  on 
the  transition  state  theory,  i.e.  a  “smoothed  potential" 
energy  surface  as  shown  in  its  diagram  of  Fig.  4  with  the 
particle  in  a  potential  well  of  barrier  height  numerically 
equal  to  the  activation  energy  of  diffusion,  Ed,  and  with  the 
width,  d,  of  the  cavity. 

The  actual  diffusion  path  along  the  diffusion  coor¬ 
dinate  proceeds  in  a  plane  normal  to  the  potential  energy 
trough  which  defines  the  saddle  area  of  width,  A  The  amount 
of  thermal  activation  Ed  required  to  allow  the  penetrant  isotope 
to  surmount  the  potential  barrier  is  (  t*  -  e^),  the  e* 
arising  from  the  partition  function  of  the  activated  complex, 

Q*,  and  from  the  partition  function,  of  the  penetrant 
in  the  cavity.  Since  ej  and  decrease  with  increasing 

mass  m^,  one  er.pects  to  observe  a  difference  in  the  measured 
isotopic  activation  energies  by  the  amount, 

6Ed  *  <Ed>D2  *<Ed>„2  "  <fD 2  *  'Da1  *  (  cH2* 

*  -<  'd2-  'h2>  +  <  {*d2-  'k2>  *  -8f  +6t*  °6) 

for  the  diffusion  of  isotopes  H2  and  D2  in  the  same  medium. 
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,  the  heat  of  solution*  is  independent  of  the  activated 
s 

state,  but  can  be  derived  from  Q^,  the  partition  function  of 
the  penetrant  in  the  cavity.  The  difference  of  the  heat  of 
solution  for  penetrants  D2  and  H2  in  the  same  medium  is 


AHs>H2  “«Da-  *H2  *  6t 


07) 


Substitute  equation  (35)  into  (3?)  for  a  mole  of  penetrant,  it 
becomes, 


Since  the  H2  and  D2  gas  particles  possess  the  same 
molecular  diameter,  0  ,  the  boundary  condition  for  the 
wave  equation  goes  in  both  cesvs  from  \<r  to  (d  -  *lo), 
yielding,  Xc, 

Replacing  by  (d— o).and  equating  d  with  one  side 
of  a  cube  which  constitutes  the  cavity,  the  incremental  free 
volume,  Vf,  one  gets 


«AHS 


-  3  N  h2 
32  (v£/3-cr  )2 


vl/3 

vf 


18,8 


+  o 


or 


(39) 
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Thus,  the  average  dimension  of  the  free  volume  can  he  estimated 
from  the  difference  in  the  heat  of  solution  for  H2  and  D£  in 
the  same  membrane,  if  the  effective  molecular  diameters  of 
the  isotopes  are  taken  to  be  equal  at  2.93  the  quantum 
mechanically  corrected  Lennard-Jones  potential  diameter  of 


H. 


(16) 


Similarily,  we  may  write 


and 


5 

P 


3  N  h2 

32 


6E_ 


3  N  h2 
32 


(ih- ;h ) 


(40) 

(41) 


since  6Ep  «6Ed  +  6AHS.  Here,  6Ep  is  the  difference  in 
the  activation  energies  for  permeation  of  two  i3otopic  species 
in  the  same  medium.  The  cavity  dimension,  Xb,  for  the  acti¬ 
vated  species  can  then  also  be  estimated. 


G.  Factors  Affecting-  P.  D  and  S 

1 .  Nature  of  the  Gas 

The  transmission  rate  of  a  penetrant  is  highly 
affected  by  the  size  of  its  molecules  since  larger  holes  need 
to  be  present  in  the  matrix  for  the  diffusion  of  larger  mole¬ 
cules  with  consequently  deeper  energy  wells.  In  fact,  in  the 
diffusion  of  bigger  molecules  lower  diffusivities  and  higher 
activation  energies  have  been  found^**®  ^  • 

2.  The  Nature  of  the  Polymer 


Because  the  diffusion  is  the  result  of  the  interplay 
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of  the  random  motion  of  the  polymer  segments  with  translational 
mobility  of  the  gas  molecules,  any  factor  influencing  the 
polymer  segmental  mobility  is  expected  to  affect  the  diffusion 

i 

coefficient.  The  presence  of  side  groups  such  as  methyls  or 
of  polar  groups,  on  the  polymer  molecule  hinders  the  chain 
movement  and  hence  should  lower  the  diffusion  coefficient  and 
increase  the  activation  energy^20  ^ •  On  the  other  hand,  the 
introduction  of  unsaturation  into  the  polymer  backbones  should 
increase  the  diffusion  coefficients  as  a  result  of  increasing 
free  volume,  or  internal  chain  mobility^1),  for  cis-conforma- 
tions  and  reduce  them  for  trans -conformations, 

3#  Effect  of  Crystallinity  or  Crosslinking 

Crosslinkings  in  the  polymer  matrix  restrict  the 

segmental  mobilities,  while  crystallites  in  the  polymer  act 

as  impermeable  obstacles  to  the  movements  of  penetrants.  Both 

effects  reduce  the  diffusivity(22»23) , 

40  Effect  of  Glass  Transition 

As  the  polymer  goes  from  the  liquid-like  rubber 

state  to  the  glass  state,  its  segmental  mobility  and  free 

volume  decrease.  The  lower  mobility  of  the  polymer  segments 

below  the  glass  transition  temperature,  Tg,  should  be  expected 

to  influence  the  permeation  process.  In  fact,  in  his  studies 

of  the  effect  of  the  glass  transition  on  the  diffusion  of 

gases  in  the  polymer,  Meares^2^  found  that  the  Arrhenius 

plot  showed  a  distinct  change  of  slope  at  T-  for  poly(vinyl 
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acetate).  However,  Stannett  and  Williams^)  found  no  inflec¬ 
tion  in  their  Arrhenius  plots  for  poly(ethyl  methacrylate). 


Kumins  and  Roteman^)  also  found  no  effect  of  the  glass 
transition  of  poly(vinyl  chloride -vinyl  acetate)  copolymer 
for  all  gases  studied  except  C02*  The  lack  of  influence  may 
be  connected  with  the  size  of  the  penetrants  and  the  relative 
extent  of  the  change  in  the  coefficients  of  isothermal  expan¬ 
sion  at  the  glass  transition  temperature.  Kumin  and  Roteman 
have  suggested  that  mainly  the  size  of  the  holes  changes  above 
Tg,  but  not  the  number.  If  the  penetrant  molecules  are  small 
enough,  the  probability  that  a  penetrant  molecule  encounters 
a  hole  of  sufficient  size  could  remain  roughly  the  same*  With 
larger  penetrant  molecules,  the  size  becomes  more  important 
and  an  effect  of  the  glass  transition  is  thus  seen^ ) , 

5.  Effect  of  Plasticizers 

The  addition  of  a  plasticizer  to  a  polymer  decreases 
the  frictional  forces  between  the  molecular  chains  resulting 
in  a  greater  segmental  mobility.  An  expected  increase  in  the 
diffusion  rate  has  been  confirmed  by  a  number  of  workers  ^"2^1 

6,  Effect  of  the  Filler 

The  permeation  rate  is  considerably  reduced  when 
the  polymer  matrix  contains  dispersed  solid  particles.  The 
effect  of  the  filler  on  diffusivity  may  be  thought  to  arise 
primarily  from  the  fact  that  each  filler  particle,  being 
impermeable,  behaves  as  a  geometric  obstruction  in  the  path 
of  the  gas  through  the  polymer.  The  gas  molecules  must  cover 
a  longer  distance  to  pass  around  the  filler  particles  and  the 
rate  of  permeation  i3  lowered*  The  shape  of  the  filler  should 
then  have  an  additional  effect  in  that  lamellar  fillers  may 
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have  a  much  greater  obstructing  geometry  than  spherical 
particles  when  the  lamellar  particles  orientate  themselves 
perpendicularly  to  the  main  path  of  the  penetrant^0 \ 

The  purpose  of  this  work  then,  as  mentioned  before, 
was  to  measure  the  transmission  rate  of  the  penetrants,  H2  or 
D2,  in  four  rubber  samples  (unfilled  and  filled  synthetic  and 
natural  rubbers)  as  a  function  of  temperatures  and  gas  pressures. 
The  transmission  rates  (or  diffusion  coefficients)  and  the 
activation  energies  for  permeation  and  diffusion,  evaluated 
from  the  Arrhenius  plot  were  to  be  compared  in  the  light  of  the 
morphological  differences  of  these  elastomers  and  with  respect 
to  the  effect  of  the  filler.  From  the  quantum  isotope  effect 
in  the  permeation :  the  cavity  size  of  the  sample  was  to  be 
estimated,  and  this  information,  coupled  with  the  densities 
of  the  samples,  used  to  evaluate  the  differences  in  physical 
structures  of  the  unfilled  synthetic  and  natural  rubbers. 
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II.  EXPERIMENTAL 


A.  Apparatus 

The  diffusion  apparatus  used  in  this  study- 
consists  of  three  main  parts \  the  gas  feeding  system,  the 
diffusion  cell  and  the  gas  measuring  system.  Figure  5  is  a 
block  diagram  for  the  apparatus. 


Fig. 5  Block  diagram  for  the  apparatus 


1.  Diffusion  Cell  and  the  Constant  Temperature  Bath 

The  diffusion  cell  was  developed  following  the  original  de¬ 
sign  of  Barrer^2),  jt  is  made  of  stainless  steel  and  consists  of 
two  parts  as  shown  in  Fig.  6,  The  upper  part  is  a  cylindri¬ 
cal  hollow  plunger  containing  a  neoprene  rubber  0-ring,  while 
the  lower  part  consists  of  a  cylindrical  hollow  base  containing 
a  metal  washer  and  a  sintered  metal  disc.  The  base  is  connected 
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Fig.  ^  Schematic  of  Diffusion  Cell 
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to  the  gas  measuring  system  via  a  Kovar  brass  to  glass 
seal.  The  plunger  is  connected  to  the  gas  feeding  system 
through  a  flexible  brass  tubing. 

In  practice,  the  specimen  is  placed  at  the  top  of 
the  sintered  metal  disc  in  the  hollow  base  which  is  in  turn 
supported  by  the  metal  washer.  Therefore,  the  specimen  will 
not  collapse  from  two  unequal  pressures  at  its  two  surfaces. 

The  O-ring  in  the  plunger  assures  a  vacuum-tight  seal  of  the 
cell  when  the  plunger  is  fastened  to  the  base  by  six  sets  of 
steel  bolts  and  screws.  The  system  is  further  leak  proofed 
by  the  mercury  seal  in  the  void  spaces  between  plunger  and 
base. 

The  constant  temperature  bath  contains  the  following 

items  i 

1)  A  fully  insulated  aluminum  bath  containing  Dow 
Corning  200  silicone  oil. 

2)  A  probe  from  the  Fisher  Temperature  Control, 
model  44,  with  a  built-in  relay  and  power  supply. 

3)  A  Bodine  motor  stirrer  from  Talboys  Instrument 
Corp# 

4)  A  quartz  immersion  heater,  lOOw. 

5)  A  Fisher  76mm  thermometer,  covering  a  temperature 
range  of  -1°C  to  110°C. 

6)  A  Fisher  thermometer  magnifier. 

7)  A  cooling  brass  coil. 

The  silicone  oil  has  e.  viscosity  of  50  cp  at  25°C.  Since  its 
viscosity  changes  only  slightly  over  a  wide  temperature  range, 
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it  should  provide  a  smooth  and  efficient  circulation  for  heat 
exchange.  By  properly  adjusting  the  speed  of  the  stirrer  and 
the  flow  rate  of  the  cooling  water*  the  bath  temperature  can 
be  controlled  to  better  than  ±0.05°C  over  the  experimental 
temperature  range,  4o-85°C. 

2#  The  Gas  veeding  System 

The  upper  part  of  the  diffusion  cell  is  connected  to 
a  two  liter  capacity  gas  storage  flask  via  a  flexible  brass 
coil  and  a  fast  action  Nupro  gate  valve  as  shown  in  Fig.  ?. 

The  flexible  coil  is  used  to  reduce  the  initial  surge  of 
gas  from  the  storage  flask  to  the  diffusion  cell  as  the  gate 
valve  is  opened.  The  storage  flask  is  connected  to  the  com¬ 
pressed  H£  and  D2  gas  cylinders  through  a  flexible  brass  coil 
and  a  stopcock.  Again,  the  brass  coil  is  used  for  reducing* 
the  surge  of  the  gas  from  the  gas  cylinder.  A  mercury  manometer 
which  records  the  storage  gas  pressure  is  at  the  left  of  the 
flask.  The  pressure  is  read  out  from  a  meter  stick  fastened 
immediately  behind  the  manometer  column.  With  the  aid  of  the 
magnifier  the  scale  can  be  read  accurately  to  ±0.1  mm  Hg.  A 
vacuum  monitor,  the  thermocouple  gauge,  is  also  incorporated 
in  the  system  so  that  the  vacuum  of  the  system  can  be  checked 
from  time  to  time.  By  means  of  a  Welch  Duo-Seal  oil  diffusion 
pump,  model  1400  (not  shown  in  the  figure),  a  vacuum  of  10“^cm 
Hg  can  be  reached  for  the  whole  gas  feeding  system.  With  the 
pump  action  isolated,  the  leak  rate  of  the  system  is  about 
one  micron  per  hour.  This  is  considered  to  be  more  than 
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sufficient  for  the  experiment,  since  the  longest  runs  lasted 
about  three  to  four  hours  and  the  shortest  runs  only  half  an 
hour. 

The  pressure  at  the  manometer' did  not  change 
during  the  course  of  experimental  runs,  since  all  runs 
were  terminated  after  the  gas  transported  had  reached  a  pres¬ 
sure  about  80  microns  in  the  receiver.  This  corresponds  to 
a  pressure  change  of  less  than  0.01$  in  the  storage  flask  for 
the  lowest  gas  pressure,  1cm  Hg»  used  in  the  experiment.  A 
manostat  was  therefore  unnecessary. 

3*  The  Gas  Measuring  System 

The  base  of  the  diffusion  cell  is  connected  to  the 
gas  measuring  system  through  a  Covar  metal  to  glass  seal  as 
shown  in  Fig.  8.  The  space  between  the  base  of  the  diffusion 
cell  to  the  stopcock  #1,  #2,  or  #3*  and  #5  determines  the 
receiver  volume  Ylf  V2  or  Vy  By  properly  manipulating  the 
stopcocks,  the  gas  transported  across  the  membrane  can  be 
routed  to  one  of  these  receiver  volumes  and  its  pressure  mea¬ 
sured  by  a  Hughes  thermocouple  vacuum  gauge  tube,  type  63^3, 
in  connection  with  ?  read-out  unit.  The  schematic  diagram  of 
the  thermocouple  gauge  is  shown  in  Fig.  9.  The  thermocouple 
(pin  3  and  5)  is  in  close  proximity  to  the  tube  heater  (pin 
1  and  7)  and  provides  an  output  B.M.F.  determined  by  the 
operating  current  to  the  tube's  heater  and  by  the  heat  conduc¬ 
tivity  of  the  residual  gas  in  the  evacuated  system.  The  ther¬ 
mocouple  output  is  10  millivolts  at  hard  vacuum,  and  0.1 
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millivolts  at  atmospheric  pressure  (air).  Thus,  this  type  of 
gauge  measures  the  change  in  output  EMF  produced  by  the  heat 
conductivity  of  the  gas  in  the  tube.  The  gas  pressure  in 
the  receiver  can  be  derived  from  the  output  EMF  through  the 
calibration  curve.  Details  of  the  calibration  procedure 
are  described  in  Sec.  II*  page  45,  An  example  of  calibra- 
tion  curve  is  shown  in  Figures  10  and  11. 

« 

The  magnitude  of  the  receiver  volumes  was  estimated 
in  the  following  manner.  A  celibration  for  the  thermocouple 
gauge  with  H2  was  established  against  the  McLeod  gauge.  The 
volume  of  the  storage  flask  in  Fig.  8  was  determined  to  be 
2003  ml  before  being  installed  into  the  system.  The  pressure 
of  H2  gas  in  and  the  volume  from  the  diffusion  cell  to 
stopcock  No.  1,  were  measured  first  and  then  a  new  pressure 
reading  was  taken  in  V2  (up  to  the  stopcock  No,  2)  with  the 
stopcock  No,  1  opened.  This  was  repeated  with  volumes, 
and  (V^  is  the  sum  of  plus  the  volume  of  the  flask) • 
The  following  volume  ratios  were  obtained  from  the  gas  lawsi 

V2  /  \1  «  1,08 

v3  /  vt  -  2.13 

v4  /  Vj  *  56.51 
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Substituting  ■  V3  +  2003  ml,  it  follows  readily  that 

Vx  »  37.9  ml 
V2  =  40.8  ml 
V3  «  80.73  ml 


to  within  approximately  sb  1%  • 

The  McLeod  vacuum  gauge,  model  2  from  Manostat 
Incorporated,  was  used  primarily  to  calibrate  the  thermocouple 
gauge.  By  assuming  that  the  gas  obeys  the  ideal  gas  law, 
pv  *  n  RT,  the  pressure  to  be  measured  is  related  to  the 
linear  length  of  the  McLeod  capillary,  L,  taken  from  the  top 
of  the  capillary  to  the  mercury  level  drawn  up  into  the  capil¬ 
lary  by, 

p  *  0.02492  L2 

where  L  is  in  units  of  mm  and  p  in  microns  of  Hg,  Since  all 
the  stopcocks  were  greased  with  Apiezon  N ,  a  high  vacuum 
grease,  and  treated  afterward  with  moderate  heat  from  an 
electrical  heat  gun  to  release  trapped  airt  the  entire  gas 
measuring  system,  including  the  base  of  the  diffusion  cell, 
was  maintained  at  a  pressure  of  1  x  10“^  cm  Kg  by  means  of 
a  mercury  diffusion  pump  and  a  Welch  DuoSeal  model  1405  fore 
pump.  The  leakage  rate  into  the  system  was  no  more  than  1 


40 


micron  ( 10 Torr)  a  day  when  the  pump  unit  was  isolated. 

B.  Materials 

1 .  Membranes 

Four  membranes  made  of  filled  and  unfilled  synthe¬ 
tic  and  natural  rubber  (SNR  and  NR),  were  used  in  this  study. 
They  were  supplied  by  the  B.  F,  Goodrich  Rubber  Go.  Akron, 
Ohio,  courtesy  of  Dr.  J.  R.  Beatty.  The  unfilled  NR  and  SNR 
(gum  stock)  have  basically  the  same  chemical  compositions, 
poly-cis-1,4  isoprene,  and  posses  the  same  low  degree  of 
crosslinking,  see  Table  1.  The  molecular  weight  between 
crosslinks  is  about  8  x  10“3#  determined  by  the  swelling  . 
method.  The  chemical  compositions  and  the  degrees  of  cross- 
linking  of  the  filled  NR  and  SNR  were  basically  the  same  as 
the  corresponding  unfilled  membranes.  The  filler  used  was 
HAF  carbon  black. 

All  the  specimens  were  cut  from  the  large  sample 
sheets,  stored  in  a  desiccator  and  kept  in  a  dark  place  at 
a  temperature  not  over  25°C  from  the  time  of  receipt.  The 
recipe  for  these  materials  is  presented  in  Table  1»  along 
with  the  pertinent  specifications  and  physical  properties. 

The  procedure  for  determining  the  bulk  density  of  the  mem-, 
branes  and  the  molecular  weight  between  the  crosslinks  is 
presented  in  the  section  of  experimental  technique  on  page  42. 


TABLE  1 
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RECIPES  AND  SPECIFICATIONS  OF  THE  SAKPIES 


Specifications  of  the  unfilled  NR  and  SNR  samples* 


Nomenclature 

(g/ cm-3) 

Crosslink¬ 
ing  agent 

Nl.wt.be  tween 
crosslinks 

Physical  Property 

Tensile  Elqnga 
strength  tion 

-  300* 

modulus 

NR  (gum  stock)** 
(p.(dsl,4-  0,9281 

isopnene)) 

dicumyl 
peroxide 
(  ~29S) 

8.28  x  103 
(by  swelling 
method  ) 

900psi  600r$ 

200psi 

SNR(gum  stock) 

(Syn.  p.  (cis  0.9183 
1,4-isoprene)) 

tl 

8.24  xlo3 

10c  170# 

•« 

♦Recipes  for  the  unfilled  NR  and  SNR  are  identical  with  the  filled 
samples  as  shown  in  the  following  Table  except  for  the  omission  of 
the  carbon  black  and  the  lower  %  dicup. 

♦♦Pale  crepe  NR 

Recipes  for  the  filled  NR  and 

SNR  samples 

(parts/weight) 

NR 

SNR 

Pale  Crepe  NR 

Dicup  40  C 
Antioxidant  2?46 
Triethanol  Amine 

N  330  Black*** 

100.0 

5.0 

1.0 

0.25 

40.0 

Syn.  p.  isoprene 

Dicup  40  C 

Antioxidant  2246 
Triethanol  Amine 

N  330  Black 

100.0 

4.25 

1.0 

0.25 

40.0 

***  N  330  black  is  an  KAF  type  (high  abrasion  furnace)  of  about 
300  A  particle  size. 
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2.  Gas»s 

Hydrogen  and  deuterium  gas  were  purchased  from 
the  Matheson  Gas  Co.  with  a  purity  of  99*95%  for  H2  and 
■>9.5#  for  D2.  The  effect  of  impurities  in  the  gases  on  the 
permeation  rates  was  negligible  and  no  further  purification 
of  the  gases  was  performed. 


C.  Experimental  Procedures  and  Technioues 


1,  Determi nation  o.f  the  Bulk  Density  of  the  Sample 
and  the  Molecular  Weight  between  Crosslinks  by 


Equilibrium  Swellin/ 


Theories  and  treatments  regarding  the  swelling  of 
the  amorphous  crossiinked  polymer  in  a  solvent  are  found  else¬ 
where  03 » 3*0 .  For  the  case  of  a  swollen  gel  in  equilibrium 

with  excess  solvent  the  molecular  weight,  j?c,  between  cross¬ 
links  is  given  by 05) 

-Pt,Vi  (sl/3-is) 

Me  - - 2 - 5 - 

In  (1-S)  +  xS*  +  s 


were  S  *  (1  +  Qg)*1#  The  ratio,  Qm(unitless) ,  of  the  volume 
of  the  liquid  in  the  swollen  gel  to  the  volume  of  the  dry 
polymer,  is  calculated  from 


pp  (Wg  -  Wp) 
PxWp 
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Subscripts,  1,  p,  and  g  refer  to  the  solvent,  dry  polymer  and 
the  gel.  Th®  densities,  p  (g/em^),  the  weights,  W,  and  the 
molar  volume  V,  were  individually  determined.  The  polymer- 
solvent  interaction  parameter,  x  ,  is  a  dimensionless  quantity 
which  is  related  to  the  energy  of  interaction  characteristic 
of  a  particular  polymer-solvent  system  and  may  be  calculated 
from  Huggin*s  equation^).  For  the  NR  (also  SNR)-benzene 
system  the  value  of  x  ■  0.425  was  taken  from  the  literature^?) 
and  used  in  the  calculation  of  Mc# 

The  procedure  used  for  swelling  the  rubber  was  the 
following^®  ^  t 

1 )  Determination  of  Wgt  To  an  amber  screwtop 
bottle  containing  O.lg  of  phenyl  naphthylamine(PNA)  and  50  ml 
of  benzene,  an  accurately  weighted  rubber  sample,  unfilled 

NR  and  SNR,  was  added  (1.5  >  0.5  gm,  approximately  dimensions 
0.3  x  0.5  x  2.5  cm).  The  bottle  was  placed  in  a  constant 
temperature  water  bath  kept  at  25  t  1°C  for  a  week.  The  sample 
was  removed  from  the  solvent,  blotted  and  rapidly  weighed  in 
a  tarred  weighing  bottle.  It  took  about  three  days  for  the 
sample  to  reach  the  swelling  equilibrium.  Duplicate  runs  were 
carried  out  for  each  sample. 

2)  Determination  of  *ulk  Density  of  the  Rubber.  Pp> 
The  bulk  density  of  the  dry  (unswollen)  rubber  was  determined 
as  follows.  Several  cubes  (each  side  about  1  mm  in  length)  of 
the  rubber  were  placed  in  a  beaker  containing  methanolc 
Water  was  then  slowly  added  from  a  burette  while  the  mixture 
was  stirred  to  insure  uniformity.  When  the  rubber  chips  were 
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suspended  within  the  mixture,  its  density  was  equal  to  that 
of  the  rubber.  The  liquid  was  then  transferred  to  a  cali¬ 
brated  pycnometer  and  weighed.  The  rubber  density  was  found 
by  subtracting  the  mass  of  the  empty  pycnometer  from  the 
total  mass,  then  dividing  by  the  volume  of  the  pycnometer. 

3)  Determination  of  the  Density  of  Solvent,  Pjj 
The  density  of  the  solvent,  benzene,  can  be  estab¬ 
lished  by  either  using  the  pycnometer  mentioned  above  or  by 
merely  talcing  literature  values  from  the  handbook.  Assuming 
a  linear  relationship  of  density  and  temperature  at  short 
temperature  intervals,  the  density  of  benzene  may  be  readily 
obtained. 

Values  of  bulk  densities  and  Mc  for  the  unfilled  NR 
and  SNR,  obtained  by  using  the  above  procedures,  are  tabulated 
in  Table  1,  in  the  material  section  on  page  41, 

2.  Sample  Preparation  for  the  Permeation  Experiment 
The  rubber  sample  as  supplied  by  the  manufacturer 
were  sheets  of  12"  x  12"  x  0.09"*  After  visual  inspection 
and  measurement  of  thickness  with  a  micrometer  this  large 
sheet  which  found  to  be  a  specimen  of  good  appearance  (no 
dent,  cracks,  scratches,  holes  or  foreign  materials)  was  cut 
in  1"  diameter.  The  thickness  of  the  specimen  was  then  mea¬ 
sured  in  at  least  ten  different  places  with  a  micrometer. 

An  average  value  of  the  thickness  was  recorded  which  was 
accurate  to  within  $  0.001  mm.  Care  was  taken  to 
avoid  surface  deformation  of  the  specimen  caused  by  exces¬ 
sive  thumbscrew  pressure  when  the  micrometer  was  applied. 
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The  specimen  was  then  mounted  into  the  diffusion 
cell  of  the  apparatus  for  degassing  in  a  manner  described  in 
Sec.  II  on  cage  30.  The  degassing  process  was  enhanced  by 
pumping  both  sides  of  the  membrane  simultaneously,  using  two 
diffusion  pump  units.  The  mercury  diffusion  pump  located  in 
the  gas  measuring  system  pumped  the  lower  side  of  the  membrane 
while  the  oil  diffusion  pump  at  the  gas  feeding  system  pumped 
the  upper  side  of  the  membrane.  This  insured  a  faster  degas¬ 
sing  rate  and  less  strain  on  the  membrane.  The  preliminary 
degassing  was  carried  out  for  two  hours  at  room  temperature. 
The  diffusion  cell  was  then  immersed  in  a  constant  temperature 
oil  bath  by  raising  the  bath  by  a  laboratory  jack.  Intense 
degassing  was  continued  at  85°C  (the  highest  experimental 
temperature)  until  the  pressure  of  the  residual  gas  escaping 
from  the  sample  was  no  greater  than  3  microns  in  half  an  hour. 
This  degassing  action  lasted  about  three  to  four  weeks  for 
unfilled  and  filled  samples.  The  need  for  this  long  degassing 
period  was  thought  to  arise  from  the  diffusion  of  residual 
solvent  and  monomers  which  were  trapped  in  the  sample  during 
the  menufacturing  process. 

3o  Calibration  of  the  Thermocouple  Gauge 

Upon  completion  of  the  degassing,  the  calibration 
of  the  thermocouple  gauge  against  the  McLeod  gauge  was  carried 
out  for  Hg  and  D2  at  room  temperature.  The  procedure  for 
calibration  was,  first,  tc  evacuate  the  low  pressure  side  of 
the  diffusion  cell,  as  well  as  the  pathway  to  the  gas  storage 
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flask,  to  the  highest  vacuum,  1  x  lCT^cm  Hg,  and  then  to  feed 
a  small  amount  of  H2  or  D2  gas  from  the  source  into  this  highly 
evacuated  space  to  contact  with  the  membrane,  thermocouple 
gauge,  and  the  McLeod  gauge.  Since  the  thermocouple  gauge  is 
relatively  insensitive  for  pressures  above  200  or  below  0.1 
micron,  the  calibration  was  not  performed  beyond  this  pressure 
range.  After  recording  the  value  of  thermocouple  gauge  output 
(in  millivolts),  the  absolute  pressure  was  determined  by  the 
McLeod  gauge.  Some  of  the  gas  was  then  allowed  to  escape 
from  the  measuring  system  and  a  new  set  of  readings  was  estab¬ 
lished.  This  was  repeated  for  approximately  ten  micron 
intervals  until  the  pressure  reached  the  final  vacuum.  This 
would  give  approximately  16  data  points  for  a  plot  of  thermo¬ 
couple  output  versus  pressure.  It  is  customary  to  establish 
the  calibration  curve  by  a  semi~Log  plot  of  the  pressure  in 
microns  of  Hg  versus  the  gauge  output  in  millivolts.  Example 
calibration  curves  are  shown  for  H2  hnd  D2  in  Figures 
10  and  11.  A  new  calibration  was  established  whenever  a 
new  membrane  was  introduced  for  an  experiment. 

4.  Collection  of  nata 

After  several  trial  runs,  the  volume  of  the  receiver, 
Vj_,  was  found  suitable  for  the  experiment.  The  runs  were  then 
carried  out  in  a  regular  order » 

1)  Work  began  with  the  highest  &s  pressure  at  the  highest 
temperature  with  one  kind  of  gas,  say  H2* 

2)  Keeping  the  temperature  unchanged  the 

gas  pressure  was  varied  from  the  highest  to  the  lowest. 


4* 

3)  The  temperature  was  then  changed  to  the  next 
lower  reading  and  the  procedure  repeated. 

4)  The  above  procedure  was  then  reputed  with 

D2  gas. 

The  degassing  time  after  each  run  depended  on  the 
temperature  and  thickness  of  the  sample.  Less  time  was  re¬ 
quired  for  thp  thinner  samples  and  the  higher  temperatures. 
For  most  experimental  condition,  about  four  hours  was  needed 
fO't*  sufficient  degassing.  However,  in  order  to  insure  a 
thorough  degassing,  more  than  four  hours  were  allowed  to 
pass  for  all  runs  before  the  data  were  collected. 

After  the  completion  of  the  degassing  stage,  "the 
permeation  experiment  was  proceeded  by 

1)  Adjusting  the  gas  pressure  in  the  storage  flask 
to  the  proper  level  by  a  bleed  valve. 

2)  Recording  the  time  with  a  stopwatch  as  soon  as 
the  Nupro  gate  valve  was  opened.  Readings  were  taken  frem 
the  thermocouple  gauge  every  half  minute  for  the  shortest  runs 
at  the  highest  pressure  and  temperature,  85°C  and  55  cm  Hg» 
and  every  five  minutes  for  the  longest  runs  at  the  lowest 
pressure  and  temperature,  1  cm.  and  40°C.  All  runs  were 
stopped  after  the  gas  transported  to  the  receiver  reached  a 
pressure  of  about  80  micron,  since  the  thermocouple  gauge 
responded  best  to  pressures  in  the  1  to  75  micron  region  (see 
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the  calibration  curve  on  pages  46  and  47 »  This  required 

about  20  minutes  for  the  shortest  runs, and  five  hours  for  the 
longer  runs, 

3)  Recording  the  gas  pressure  accurately  from  the 
monometar  with  the  aid  of  the  4"  magnifier,  and  isolating  the 
gas  source  from  the  diffusion  cell  by  closing  the  Nupro  gate 
valve, 

4)  Finally,  degassing  the  membrane  for  at  least 
four  hours  before  the  next  run  was  performed, 

5*  Calculation  of  P,  D,  and  S 

If  the  gas  obeys  the  ideal  gas  law,  pv  *  n  RT,  the 
pressure  in  the  receiver,  Vj_,  may  be  treated  as  proportional 
to  the  amount,  Q^,  transported.  Applying  equation  (15),  a 
plot  of  Q-c  versus  t,  was  made  for  each  of  the  permeation  runs. 
The  time  lag,  t  was  obtained  by  extrapolating  the  linear 
portion  of  the  plot  to  the  time  axis  as  shown  in  Figure  12 
while  che  slope  of  the  line  gave  the  flux,  Qs,  However, 
this  graphical  extrapolation  to  the  time  axis  is  equivalent 
to  evaluating  the  zero  time  for  the  linear  pressure  rise 
from  the  ideal  zero  value.  In  order  to  take  into  account 
the  small  residual  pressure  in  the  receiver,  V1#  at  the  be¬ 
ginning  of  each  run,  the  following  analytical  method  was 
used  to  obtain  the  value  of  r  rather  than  merely  reading 
the  intercept  at  the  time  axis. 

The  linear  portion  of  the  curve  in  Fig.  12  can 


be  represented  by 


3K5iIXCTl5iWF.ffi3rcfc  £ 
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p  *  a  +  bt 


where  p  is  the  pressure  in  microns,  b  is  the  .slope  in  micron 
per  minute,  t  is  the  time  in  minutes,  and  "a”  is  the  intercept 
in  microns.  Solving  the  equation  for  t  *  0,  at  p  *  p0,  one 
obtains 


min  - 


The  diffusion  coefficient,  D,  is  calculated  from  the  value  of 
r  as  shov/n  in  equation  (16) 


l/x  b 


cm2/s 


The  slope  of  the  line  in  Fig.  .12  represents  the  steady  state 
flow,  Qs,  through  a  cross  sectional  area,  A.  The  permeability, 
P»  in.  equation  (9)  is  then  given  by 


Os  x  L 

p  ■  - - 

A  *  p 


(43) 


The  solubility,  S,  follows  readily  according  to  fcjuatior.,  (11), 
S  «P/D.  A  sample  calculation  of  P,  D  and  S  is  given  in. 
Appendix  I , 


6.  Calculation  of  Ed,  E^, 


and  AHs 

jThe  Arrhenius  equations  (25),  quoted  on  page  16,  maj, 
be  written  in  the  logarithmic  form  as 

log  D ■ log  D°  -  ?  (44) 
and 

log  P  =  log  PQ  -  2.303Pr  «jf  (^5) 

where  the  value  of  R  is  1.98?  cal/mole -degree.  Thus,  a  plot 
of  log  D  or  log  P  against  l/T  should  give  a  straight  line  ,  the 
slope  of  which  equals  -Ed/2.303R  or  -Ep/2.303R  .  The  activa¬ 
tion  energy,  Ed,  or  Ep,  can  thus  be  evaluated. 

However,  in  order  to  obtain  a  best  straight  line,  the 
graphical  method  was  not  used.  Instead,  the  least  square  irathod 
was  employed  to  calculate  the  slope  and  the  intercept  of  a  best 
straight  line  from  the  data.  Since  both  equations  above  can 
be  simply  represented  by 

y  *  a  +  bx  (46) 

where  y  *  log  D  or  log  P 

a  =  log  D0  or  log  P0 
b  *  Ed/2.303R  orEp/2.303  R 
x  =  l/T 

A  sample  calculation  of  Ed  and  Ep  is  shown  in  Appendix  II, 

AHS  follows  as  the  difference i  (Ep  -  Ed) . 


Ill  EXPERIMENTAL  RESULTS 


This  section  is  divided  into  four  parts,  devoted 
to  the  fo’T  types  of  samples.  Values  of  permeability,  P, 
diffusivity,  D  and  solubility,  S  at  various  gas  pressures 
and  temperatures  studies  for  each  sample  are  calculated  in 
the  same  manner  as  shown  in  the  '‘sample  calculation"  in 
Appendix  I  and  tabulated  accordingly.  The  Tables  of  P,  D, 
and  S  are  proceeded  by  a  summary  of  the  activation  energies 
for  permeation  and  diffusion,  Ep  and  Ed,  the  heat  of  solution, 
AHs,  and  the  pre -exponential,  D0*  at  each  gas  pressure.  Ed, 

Ep  and  D0  are  calculated  by  the  least  square  method.  A  sample 
calculation  i*a  given  in  Appendix  IT.  AHS  is  the  difference 
of  Ep  and  Ed. 

In  all  cases,  the  area  of  the  membrane  exposed  to 
the  penetrant  was  2.35  cm2  cell  and  a  receiver  volume,  Vj_,  of 
37.9  ml  was  used.  The  average  thickness  of  each  sample  is  as 
the  follows  t 


Sample 
Unfilled  NR 
Unfilled  SNR 
Filled  NR 
Filled  SNK 


Average  thickness,  cm 
0.2101 
0.2116 


0.2097 

0.1926 


'l  i 


The  Arrhenius  plot  is  practically  linear  for  the  temperature  range 
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IV  DISCUSSION  OF  RESULTS 


The  experimental  results  listed  in  Table  2  through  13 
in  Section  III  are  discussed  in  the  order  of  (A),  diffusivity, 
D,  solubility,  S,  and  permeability,  Pt  (B),  activation  energies 
for  diffusion  and  permeation,  Ed  and  Ep,  and  heat  of  solution, 
£HS  and,  (C),  the  free  volume  of  rubber.  Except  for  the  free 
volumes,  all  the  above  quantities  are  treated  as  a  function  of 
gas  pressures  and  in  one  case  as  a  function  of  temperature. 

A.  Diffusivity,  D,  Solubility,  S,  and  Permeability,  P, 

1.  Diffusivity,  D,  as  a  Function  of  Gas  Pressure 

The  diffusion  coefficients  evaluated  from  the  time 
lag  r,  according  to  equation  (16),  assuming  a  constant  D, 


T 


for  the  transport  of  hydrogen  and  deuterium  gas  through  the 
four  samples,  are  separately  plotted  against  the  gas  pressure 
as  shown  in  Figures  13  to  16,  In  examining  these  figures,  one 
finds  that  these  diffusivity  isotherms  are  indeed  largely  inde¬ 
pendent  of  the  gas  pressure  at  low  temperatures,  and  become 
slightly  pressure  dependent  at  higher  temperatures  but,  in  no 
case,  does  the  difference  of  the  diffusivities  through  the  pres¬ 
sure  range,  i.e.,  between  1  and  55  cm  of  Hg,  for  any  isotherm 
exceed  the  unit  value  (  lcm2/sec  )  of  diffusivity.  The  average 
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values  of  each  isotherm  are  listed  in  Table  14* 

Comparison  of  these  values  reveals  three  effects  which 
influence  the  magnitude  of  the  diffusivlty  (besides  the  tempera¬ 
ture)!  the  effects  of  the  penetrant  mass*  that  of  the  nature 
of  the  rubber i  and  the  effect  of  the  filler  in  the  rubber*  The 
first  indicates  that*  in  the  same  medium*  hydrogen  has  a  higher 
diffusivity  than  its  isotope*  deuterium t  the  second  shows  that 
this  difference  varies  very  much  with  the  nature  of  the  rubber! 
and  the  third  effect  indicates  that*  for  the  same  penetrant* 
the  diffusivity  in  the  filled  rubber  is  several  times  lower 
than  in  the  unfilled  rubber* 

(1)  Effect  of  the  penetrant  massi  An  influence  of 
mass  is  obvious  if  one  considers  the  translational  kinetic 
energy  of  a  penetrant  along  the  diffusion  coordinate,  yielding 

Efc  »  J  kT  ■  £  raivi^ 

This  equation  shows  clearly  that  at  the  same  temperature*  the 
lighter  penetrant  should  have  a  higher  penetrating  velocity* 

For  Hjj  and  D2  the  mass  ratio  is  2,  and  therefore  the  diffusion 
ratio  should  be  1*4*  For  unfilled*  and  less  so  for  filled*  NR 
we  see  the  ratios  to  be  close  to  the  ideal  values*  within  the 
experimental  error*  For  SNR*  there  is  practically  no  difference 
between  the  diffusivities  of  H2  end  D2*  This  is  quite  remarkable 
because  the  difference  in  mass  must  affect  the  diffusion  mechanism. 
There  must  be  a  compensating  mechanism  that  enhances  the  diffusion 
of  D2  in  SNR,  possibly  connected  with  the  lack  of  crystallinity  in 
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TABLE  14 

AVERAGE  VALUE  OF  DIFFUSIVITIES  FOR  GAS  PRESSURES  IN  THE 
RANGE  OF  1  TO  55  cm  Hg  AT  CONSTANT  TEMPERATURE 


Gas 

t 

Unfilled  NR 

Filled  NR 

Ratio 

<  °C) 

D 

D 

(°unf./  Dfilled) 

40 

1.63 

0.39 

4.19 

55 

2.45 

0.65 

3.78 

H2 

70 

3.24 

1.02 

3.19 

85 

4.57 

1.61 

2.85 

40 

1.18 

0.35 

307 

55 

1.76 

0.58 

3.06 

d2 

70 

2.51 

0.88 

2.84 

85 

3.44 

1.29 

2.66 

t 

Unfilled  SNR 

Filled  SNR 

Ratio 

Gas 

(  °C) 

D 

D 

(°unf ./  ^filled) 

40 

1.71 

0.45 

3.78 

6o 

2.82 

0.86 

3.27 

h2 

70 

3.51 

1.12 

3.12 

85 

4.58 

1.72 

2.66 

40 

i..6o 

0.43 

3.75 

50 

2.18 

0.57 

3.82 

d2 

60 

2.8? 

0.76 

3.77 

70 

3.35 

0.99 

3.37 

85 

4.63 

1.44 

3.21 

D  in 

cm2/sec. 

is  0.1 6  volume 

Filler 

content 

fraction. 
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the  latter,  £ee  also  Table  17  )• 

(2)  Effect  of  filleri  Fillers  contained  in  the  rubber 
may  be  reinforcing  or  non-reinforcing  and  their  particles  may 
be  spherical  or  non-spherical  depending  on  the  type.  The  most 
important  groups  are  mineral  fillers,  fibrillar  or  lamellar 
fillers,  reinforcf  ~  and  non-reinforcing  carbon  blacks,  and  organic 
fillers  such  as  resins,  pigments  or  other  rubber  phases.  The 
particular  type  of  filler  (carbon  black)  used  here  is  a  reinfor¬ 
cing  one.  Consequently,  the  filler  particles  will  not  only 
act  as  impermeable  barriers  but  also  affect  the  structure  of 
the  rubber  in  their  vicinity  b,  various  types  of  rubber-filler 
attachments.  If  the  surfaces  of  the  particles  are  wetted  by  the  rubber 
and  thus  induce  a  more  or  less  dense  distribution  of  rubber  chain 
segments  at  the  interface,  the  filler  particles  will  certainly 
have  the  effect  of  reducing  the  overall  chain  flexibility.  This 
is  evid-anced  bydatain Table  14.  The  ratio  of  the  average  diffusivity 
for  the  unfilled  rubber  to  the  filled  ranges  approximately  from 
2.6  to  4,2,  and  the  ratio  decreases  as  the  temperature  increases. 
The  latter  fact  implies  that  the  diffusivity  increases  faster 
with  the  temperature  in  the  filled  rubber  than  in  the  correspond¬ 
ing  unfilled  one,  probably  as  the  result  of  a  desorption  of  the 
locally  concentrated  segments  from  the  filler  surfaces,  thus 
making  the  filled  rubber  relatively  more  flexible. 

Diffusivities  measured  by  van  Amerongen^O  *31)  for 

natural  rubber  containing  various  types  of  fillers  are  quoted 

(40) 

in  Table  15.  Also  of  interest  are  the  data  of  Ziegel  • 
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Comparing  the  uncorrected  value  for  D  of  the  unfilled  samples  , 
without  adjusting  for  the  difference  in  temperatures ,  we  find 
D  =  1.63  x  1CT5  versus  D  =  1.85  x  10"^  from  von  Amerongen's 
data.  This  is  in  excellent  agreement  if  one  takes  into 
account  the  three  degrees  difference  in  temperature,  and  the  in¬ 
evitable  differences  in  rubber  samples.  Taking  the  filled  sample?, 
a  comparison  becomes  more  difficult.  But  we  can  see  that,  for 
the  comparable  carbon  black,  Amerongen  finds  the  ratio  of  the 
unfilled  to  the  filled  samples  to  be  also  in  the  range  of  4  to 
5.  Comparing  with  Ziegel,  his  value  of  1.82  x  10“^  at  44°C  for 
the  unfilled  SNR  for  deuterium  is  also  in  excellent  agreement 
with  our  value  of  1.60  x  10~5  at  40°C.  For  hydrogen  there  is 
also  good  agreement.  As  shown  by  Meares^^  ,  these  values  of 
diffusivities  are  typical  for *  hydrocarbon  rubbers.  The  present 
values  for  the  filler  effect  are  also  supported  by  other  relevant 
data  in  the  literture^2'^) , 

The  present  data  differ,  however,  from  Ziegel's  in 
that  they  do  not  show  the  strong  pressure  dependence  of  the 
diffusion  constant  at  low  pressures.  It  is  likely  that  in  this 
range  Ziegel's  data  exhibit  a  spurious  trend  because  his  thermo¬ 
couple  gauge  was  not  as  accurate  as  the  one  used  here#  From  the 
gas  pressure  of  about  10  cm  Hg  on,  our  and  his  data  are  quite 
comparable.  ' 

Concerning  the  magnitude  of  the  filler  effect,  we  see 
that  permeation  has  been  reduced  much  beyond  the  change  in  pene¬ 
trable  cross-section  which  is  about  36%  of  the  unfilled  rubber. 

In  fact,  tne  change  in  D  is  greater  than  the  factor  of  about  2 
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by  which  our  volume  fraction  of  filler  would  increase  the  rubber 
viscosity,  or  modulus,  (assuming  Du  *  constant)  in  accordance 
with  any  of  the  existing  equations  for  this  effect,  such  as 
Mooney's^^),  Guth's^^,  Smallwood’s^^,  or  Kerner's^^  equation. 
In  other  words,  the  reduction  in  diffusivity  must  be  due  tc  yet 
another  factor  besides  the  corresponding  reduction  in  the 
mobility  of  the  rubber  between  the  filler  particles. 

In  the  absence  of  a  detailed  theory  of  the  change  in  diffusion 
due  to  the  replacement  of  rubber  by  barrier  particles,  it  would 
be  difficult  to  estimate  how  much  of  this  is  due  to  the  changes 
in  bulk  rubber  mobility  and  how  much  due  to  the  strain  or  com¬ 
paction  of  the  rubber  adsorbed  to  the  filler  particles.  The  tor- 
tuosityaf  the  diffusion  path  is  also  involved  .  The  situation  is  actually 
similar  to  the  diffusion  through  a  partially  crystalline  polymer 
where  the  crystalline  volume  fraction  acts  like  a  filler.  This 
will  be  discussed  in  later  Sections. 

2.  Diffusivity  as  a  Function  of  Temperature 

The  isobaric  plots  of  diffusion  coefficients  for  our 
four  samples  are  shown  ;n  Figures  1?  through  20.  In  all  the 
plots  the  degree  of  divergence  indicates  the  pressure  dependence 
of  D  which  again  is  seen  to  be  small  at  low  temperature,  but 
increases  at  higher  temperature. 

The  plots  for  the  unfilled  samples  are  linear  in 
contrast  to  those  of  the  filled  rubbers.  The  linear  isobaric 
plots  (for  the  unfilled  rubber)  may  be  represented  empirically 
by  the  following  four  equations! 
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D2-unfilied  SNR  D  =  6. 19  (t2-t1)  x  10“2  14.5  t  5,7% 


D  in  unit  of  cm2/sec  x  10“5 


Here,  denotes  a  practical  "cut-off"  temperature  at  which  the 
diffusivity  appears  to  vanish  (in  fact,  becomes  very  small  com¬ 
pared  with  values  between  40  and  85°C).  Using  this  linearized 
equation,  the  diffusivity  at  any  temperature,  t2»  can  be  calcu¬ 
lated  in  this  range  with  an  error  of  no  more  than  6,8#  compared 
with  the  average  value  listed  in  Table  14.  However,  diffusivi- 
ties  calculated  for  temperature  higher  than  the  experimental 
upper  bound  temperature  of  85°C  may  be  sub ject  to  greater  errors. 

The  "linearity"  of  D  with  temperature  in  the  unfilled 
samples  expresses  the  fact  that  the  diffusion  mechanism  resembles 
that  of  a  diffusant  following  Einstein’s  law  if  moving  through 
a  medium  whose  viscosity  changes  little  with  the  temperatures 
D  3  kT/6irnr.  The  real  curvature  of  those  plots  is  not  only 
implied  by  the  applicability  of  Arrhenius  plots  which  will  be 
discussed  later,  but  becomes  strongly  pronounced  in  the  filled 
samples,  where  rubber  attachments  and  viscosities  must  change 
more  strongly  with  temperature. 
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3*  Solubility.  S,  as  a  Function  of  Gas  Pressures 

The  solubilities  listed  in  Tables  3*  4-»  6,  ?,  9»  10,  12 
and  13  are  plotted  against  the  gas  pressure  for  two  temperatures, 
40  and  85°C  (the  lower  and  upper  bound  temperature  used  in  the 
study)  as  shown  in  Figures  21  to  24.  It  is  interesting  to  note 
that  in  both  unfilled  rubbers,  the  two  penetrants  show  the  same 
trend  of  increasing  solubility  with  temperature  whereas  in  the 
filled  rubbers  the  two  isotope  species  show  different  trends  in 
their  solubilities.  For  hydrogen  the  solubility  in  the  filled 
rubbers  decreases  as  temperature  increases,  while  for  deuterium 
the  trend  is  reversed  and  less  regular. 

These  differences  can  be  explained  by  the  signs  of 
their  heats  of  solution.  From  the  thermodynamic  point  of  view, 
positive  values  of  heats  of  solution,  i,e.„  endothermic  dissolu¬ 
tion  of  the  gas,  cause  rising  solubilities  with  increasing 
temperatures.  On  the  other  hand,  for  negative  heats  of  solution, 
or  exothermic  processes,  the  solubility  should  decrease  as 
temperature  increases.  The  mean  heats  of  solution  deduced  for 
the  respective  penetrant-rubber  systems  listed  in  the  following 
(quoted  from  Table  l6  on  page  125  )  are  consistent  with  this, 
as  shown  be low t 


^Hof  kcal/mole 

o 


^v.hubber 

Unfilled 

Filled 

NR 

SNR 

NR 

SNR 

h2 

1.24 

1.02 

-0.28 

-0.34 

D2 

0.65 

0.71 

0.28 

0.35 

GAS  PRESSURE  ,  CM  OF  HG 


23  -  Solubility  isotherms  versus  gas  pressure 
for  H2  and  synthetic  natural  rubber 
(filled  or  unfilled). 
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D  -  Filled  SNR 


4^40°C 


D2  -  Unfilled  SNR 

_a _ A 


5^40  50  60 


10  20 
GAS  PRESSURE  .  CM  OF  HG 

«• » -  r  sf  s  «S“r““’”s“r 

(filled  or  unfiUe0>* 


In  this  Tabic,  AH  is  the  difference  of  the  arithmetic  mean  of  the 
activation  energies  for  permeation  and  for  diffusion  for  the 
pressure  range  investigated  (see  Section  of  activation  energies 
on  page  117  ).  The  value  of  AHg  for  each  gas  pressure  can 
also  be  calculated  from 


*Hs 


-R  d  InS 
(1/T) 


which  is  a  more  direct  method.  However,  if  this  is  done,  the 
trends  of  AHS  are  slightly  clearer  but  the  scatter  of  the  values 
remains.  Therefore,  values  of  AHS  were  calculated  by  the 
former  method. 

The  solubility  isotherms  shown  in  these  figures  also 
illustrate  the  fact  that  the  solubilities  of  the  penetrants  in 
the  filled  rubbers  are  much  higher  than  in  the  unfilled  rubbers. 
In  general,  the  solubility,  s(rp+f)»  of>  a  6as  in  a  rubber-filler 
system  may  be  resolved  into  three  factors,  first  Srb’  related  to 
the  solubility  in  the  rubber,  second,  Sf,  related  to  the  gas 
adsorption  by  the  f.iiier  and,  third,  5vac,  related  to  filling 
up  (condensing  in)  vacuoles  or  small  capillaries  which  may  be 
present  in  the  porous  fillers,  or  may  be  created  by  high  concen¬ 
tration  of  fillers  in  the  rubber^  58) 

s(rb+f)  =  xrbsrb  +  xfsf  +  XvacSvac 

where  Xrt)#  Xj*,  and  *vac  are  the  volume  fractions  of  rubber, 
filler  and  vacuolts  respectively.  In  the  special  case  of  com¬ 
plete  wetting  of  the  surface  of  a  non-purous  filler  by  the  rubber, 
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the  above  equation  reduces  to 

s(rb+f)  3  xrbsrb 

With  rubber  containing  filler  of  say,  0.2  volume  fraction,  the 
solubility  of  the  gas  should,  in  this  case,  be  80#  of  that  of 
the  gum  compound.  For  special  mineral  fillers  such  as  whiting, 
aluminum  oxide  and  barium  sulfate  and  with  the  non-reinforcing 
carbon  blacks,  the  experimentally  determined  solubilities  by 
van  Amerongen^0)  are  in  fact  only  a  little  higher  than  this 
80#  value  (see  Table  15)#  If  the  filler  particles  contained  in 
our  filled  rubber  were  non-adsorbing,  the  solubilities  of  the 
gas  should  be  about  84#  of  that  of  the  unfilled  rubber  since 
the  volume  fraction  of  the  filler  is  0.16.  In  fact,  we  find 
for  our  system  that  is  almost  twice  as  much  as  Srb# 

This  is  certainly  due  to  the  high  adsorption  power  of  the 
carbon  black  and  also  possibly  due  to  the  increase  in  free  volume 
as  a  result  of  the  strain  introduced  by  the  filler  particles. 

The  later  effect,  of  course,  should  be  small. 

Figures  21  to  24  show  interesting  trends  of  solubili¬ 
ties.  The  solutions  of  H2  gas  in  the  filled  rubbers  are  the 
only  ones  showing  the  negative  AHg  values.  These  are  undoubtedly 
due  to  the  strong  adsorption  of  hydrogen  on  the  filler  and  thereby 
presents  special  cases.  However,  even  these  two  (like  the  other 
six  cases)  show  a  general  trend  of  increase  in  solubilities  at 
low  pressures. 
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One  sees  that  S  is  reasonably  constant#  i.e., 
that  Henry’s  law  is  fulfilled,  for  three  of  the  unfilled 
rubber  experiments.  In  the  forth  case,  there  is  a  downward 
trend  of  solubility  with  pressure  which  is  probably  exag¬ 
gerated  by  experimental  errors.  In  the  filled  rubbers, 
the  solubility  declines  quite  markedly  with  pressure,  par¬ 
ticularly,  at  the  low  pressures,  which  indicates  compli¬ 
cations  beyond  the  ordinary  solubility  processes.  This  may 
derive  from  a  superimposed  effect,  according  to  Gibbs’  equa¬ 
tion,  if  a  gas  acts  as  external  pressure  carrier  besides 
dissolving  in  the  material  under  pressure.  Very  likely 
we  see  a  variation  of  this  in  the  form  of  an  effect  of  me¬ 
chanical  strain,  intensified ‘by  the  presence  of  fillers  and 
vacuoles,  as  the  rubber  is  compressed  against  the  porous 
net&l  support  in  the  diffusion  cell. 

The  data  found  here  for  the  solubility  and  the 
heat  of  adsorption  allow  an  interesting  calculation.  Assuming 
an  average  surface  area  of  the  carbon  (particle  size  300  A  in 
disaster)  of  120  m  /g,  and  9  A  for  the  area  required  for 
each  H2  nolecule,  one  can  calculate  that  the  amount  of  the 
carbon  required  for  adsorption  of  one  mole  of  H2  gas  is  about 
500  g.  In  view  of  the  weight  ratio  of  the  rubber  to  the  carbon 
in  the  filled  rubber,  there  should  be  about  2.5  times  as  much 
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rubber  required  to  adsorb  one  mole  of  H£  gas.  Instead,  our 
solubility  data  at  gas  pressure  of  30.2  cm  Kg  and  40°C  reveal 
that  only  8  x 10*^  g.  of  H2  gas  is  adsorbed  on  1  g.  of  carbon 
in  the  filled  NR.  In  other  words,  2.5x10^  g,  of  carbon  is 
needed  to  adsorb  one  mole  of  H2  gas  in  our  case.  The  ad¬ 
sorptive  capacity  of  the  carbon  in  the  filled  rubber  is  thus 
only  1/500  of  that  calculated  for  a  theoretical  monolayer. 
Undoubtedly,  this  is  a  result  of  a  lower  effective  surface 
area  of  the  carbon  due  to  the  fact  that  probably  90 ft  or 
more  of  the  carbon  surface  is  occupied  by  the  rubber. 

The  experimental  heat  of  adsorption  for  H2  on  carbon 
is  -1.52  kcal/mole  which  is  obtained  by  taking  the  difference 
of  AHS  for  filled  and  for  unfilled  NR.  The  heat  of  internal 
adsorption  for  H2  on  rubber  is  estimated  about  -0.8  kcal/mole. 
This  is  reasonable  since  this  energy  must  be  smaller  than 
the  heat  of  adsorption  on  carbon  and  larger  than  the  product 
of  entropy  loss  and  temperature  due  to  the  gas  being  adsorbed 
at  low  pressure  and  room  temperature.  This  loss  is  about 
600  cal/mole,  i.e.,  RT.  Therefore,  the  energy  of  hole 
formation  is  found  to  be  2,040  cal/mole  (see  the  analysis  of 
the  heat  of  solution  on  page  126,  where  the  heats  of  self 
condensation  of  gas  in  the  rubber  and  on  the  carbon  are 
stated  to  be  negligibly  small).  This  is  about  two  fifth  of 
the  frequently  derived  molar  cohesive  energy  of  rubber  and 
again  reasonable. 

The  following  consideration  is  also  interesting^  The 
heat  of  solution,  AHS  for  H2  in  unfilled  NR  is  1.24  kcal/mole. 
To  occur  spontaneously,  the  free  energy  of  the  solution  must 


be  negative.  According  to  Gibbs*  equation,  AF  *  AHS  -  TAS , 
the  entropy  ofmixing,  AS,  must  be  positive  and  the  term,  TAS, 
must  be  greater  than  AHs  which  is  1,24  kcal/mole,  in  order  to 
have  a  negative  AF.  Again  setting  the  entropy  loss  due  to  H2 
adsorbed  on  rubber  to  be  -2.0  e.u/mole  and  T  *=  300°K,  there 
must  be  an  entropy  increase  in  the  rubber  of  at  least  6  e.u 
per  mole  of  isoprene  through  H2  adsorption.  This  "entropy  of 
swelling"  is  large,  but  at  least  comparable  to  the  entropy  of 
fusion  of  NR,  3*5  e.u/mole,  quoted  in  the  literature ^ . 

4.  Permeability,  P 

The  values  for  P  found- have  been  given  in  Tables  3» 
4,  6,  7.  9,  10,  12,  and  13.  Obviously,  the  magnitudes  and 
trends  of  permeabilities  agree  with  the  discussion  above  on 
the  factors  that  influence  D  and  S,  of  which  P  is  the  product. 
However,  one  must  remember  that  the  permeabilities,  as  norma¬ 
lized  flux  rates,  are  the  primary  data  since  even  the  values 
of  D  are  derived  from  the  intercept  of  the  flux  rate  on  the 
time  axis.  Consequently,  features  of  the  membrane  structure, 
or  of  the  absorption  behaviour  of  the  gases,  must  receive 
here  some  attention. 

Concerning  the  unfilled  rubbers,  Ziegel^0  ^  e.g*, 

noted  a  change  in  permeability  of  the  samples  with  time.  He 
attributed  this,  probably  correctly,  to  aging  and  crystalliza 
tion  of  the  rubber  membranes,  facilitated  by  the  changes  in 
temperature  and  pressure,  and  the  presence  of  the  dissolved 
gas.  Such  time  effects  on  permeabilities  were  not  observed  in 
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this  study. 

On  the  other  hand,  in  dealing  with  the  filled  samples 

one  must  analyze  the  possible  changes  of  rubber  structure  and 

diffusion  pattern  due  to  the  presence  of  the  reinforcing  carbon 

black.  Firstly,  there  is  a  change  in  the  effective  cross 

section  for  the  diffusive  passage  of  the  gas  v/hich  the  filler 

reduces  to  8 kfi  of  the  total  membrane  area.  Secondly,  we  have 

to  consider  the  fact  that  the  gas  molecule  can  not  follow  a 

direct  diffusion  path,  but  must  bypass  the  filler  particles, 

so  that  the  effective  travelling  length  through  the  membrane 

becomes  greater.  This  so-called  tortuosity  factor  has  been 

repeatly  discussed  in  the  literature  in  connection  with  the 

flow  cf  fluids  through  porous  media.  Although  the  permeation 

of  gas  through  the  membranes  is  a  matter  cf  diffusive  rather 

* 

than  fluid  mass  transport,  the  generally  accepted  tortuosity 
factor  of  1,6  for  randomly  packed  beds^)  should  apply,  though 
as  a  maximum.  Thus,  the  effective  membrane  thickness  as  greater 
by  factor  of  1.6,  or  the  normal  concentration  gradient  estab¬ 
lished  in  the  membrane  should  be  reduced  by  a  factor  of  1.6. 

Thirdly,  the  mobility  of  the  rubber  segments  located 
next  to  the  filler  particles  must  be  reduced,  bringing  the 
segmental  viscosity  closer  to  that  prevailing  in  the  glassy 
state.  The  observation  that  temperature  has  such  a  large 
effect  on  the  diffusivity  of  the  filled  rubbers  could  be  ex¬ 
plained  on  this  basis. 

Finally,  the  increased  solubility  of  the  gases  in 
the  filled  samples,  as  well  as  the  overall  more  negative 


enthalpy  of  solution,  suggests  higher  concentrations  of  the 
gases  in  and  around  the  filler  particles.  Consequently,  there 
may  be  another  mechanism  of  diffusion  besides  the  assumed 
migration  of  the  gas  molecules  from  segmental  site  to  site  by 
adsorption  and  desorption  steps  on  the  rubber  This  latter 
mechanism  neglects  another  likely  aspect,  namely  that  rubber 
segments  themselves  also  diffuse.  In  doing  so  they  enhance  gas 
transport  by  moving  out  of  the  way  of  gas  particles,  thereby 
opening  suitable  cavities  for  the  gas  to  diffuse  into,  at  an 
energy  expense  which  is  different  from  that  of  the  activated  desorption 
assumed  in  this  thesis.  This  new  mechanism  would  allow  the 
gas  molecules  to  diffuse  into  the  dense  gas  zone  in  and  around 
the  filler  particles,  to  become  adsorbed  there,  while  a  gas 
particle  at  the  opposite  side  of  the  filler  is  given  off. 

This  involves  a  fluid  type  of  transport  in  the  adsorbed  atmos¬ 
phere  of  the  filler  which,  if  this  movement  is  effective, 
would  decrease  the  obstructive  power  of  the  filler.  There  may 
be  yet  other  factors  at  work,  but  we  want  to  restrict  the 
discussion  to  these  obvious  ones. 

The  combined  effects  of  change  in  effective  membrane 
cross  section  through  the  filler  and  of  rubber  adsorption  on 
the  carbon  may  be  described  as  follows i  We  accepted  the  rela¬ 
tion  P  *  SD,  as  a  result  of  Fick's  first  law,  0  *  *  and 

of  Henry’s  law,  c  =  Sp,  under  the  assumption  of  constant  D  (see 
derivation  on  page  6  ).  However,  since  the  cross  section  is 

changed  from  unity  to  an  area  reduced  by  to  8^,  and  since  the 
length  of  the  diffusion  path  is  increased  by  6 0$,  one  sees 
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that  the  permeability,  P,  for  the  filled  rubbers  must  be  re¬ 
duced  to  one  half  of  that  for  the  unfilled  rubber  for  these 
geometrical  reasons  alone.  Introducing,  then,  the  effect  of 
fillers  on  the  segmental  mobility  arid  the  consequently  increased 
viscosity,  as  discussed  in  the  section  of  diffusivity,  we  first 
arrive  at  a  diffusion  rate  ratio  between  the  unfilled  and 
filled  samples  of  0,55*  However,  the  time  lag  itself  must 
also  be  affected  by  tortuosity  so  that  the  diffusion  factor 

p 

would  be  further  reduced  by  1,6,  or  by  (1.6)  as  maintained 
>  f  *0 ' 

by  some  theories'-^  ' .  Applying  these  factors,  leads  to  the' 
diffusion  ratios  as  listed  in  'rable  14. 

It  follows  from  this  analysis  that  our  values  for 
solubilities,  S,  in  the  filled  samples  are  also  affec-ted.  The 
true  solubility  should  be  the  quotient  of  the  true  permeability 
divided  by  the  true  diffusivity.  The  true  permeability  is  in 
reality  larger  than  that  measured  if  one  takes  the  tortuosity 
and  the  change  in  the  cross  section  of  the  specimen  into 
account.  Similarly,  the  true  diffusivity  should  be  cor¬ 
rected  for  the  tortuosity  (  or  the  true  length  of 
the  diffusion  path  )•  If  all  these  corrections  are 

made,  the  true  solubility  calculated  would  be  smaller  than  the 
apparent  value  (i.e*,  calculated  from  the  uncorrected  values 
of  P  and  D)  by  a  factor  of  1.35*  No  doubt,  this  argument  can 
be  refined  or  revised  by  a  different  model  of  filler  action. 

For  the  present,  we  have  calculated  the  solubility  directly 
from  the  experimental  permeabilities  and  diffusivities  which 
are  uncorrected  for  the  tortuosity.  Having  thus  preserved 


the  primary  data  they  can  be  later  recalculated  in  accordance 
with  a  more  complete  theoretical  study  of  filler  action. 
Similarly,  the  quantities  of  Hz  adsorbed  on  carbon  which  we 
have  derived  with  the  aid  of  solubilities  might  have  to  be 
revised  later.  But  even  without  this  correction,  the  essence 
of  the  conclusions  drawn  here  remains  correct. 

The  temperature  dependencies  of  the  permeabilities 
and  diffusivities  can  be  expressed  by  the  Arrhenius  equation. 
Values  at  different  temperatures,  T2,can  be  calculated  from  a 

value  of  Ti  and  the  activation  energies  by  the  following 
relationships! 

log  P2  =1=6  pt  ♦&(£-£) 

Ed  /  1  1  \ 

log  D2  -  log  D2  +  — (  -  -  -  ) 

where  Pj ,  F2  and  Dj.,  D2  are  the  permeabilities  and  diffusivi¬ 
ties  at  temperatures  T*  and  T2  °K  respectively,  and  Ed  and  Ep 
are  the  activation  energies  in  calories  per  mole,  see  Figures 
25  and  26,  which  are  presented  as  examples. 

The  temperature  dependencies  of  the  diffusivity  and 
permeability  are  usually  such  that  P  and  D  decrease  uniformly 
with  decreasing  temperatures.  However,  in  some  temperature- 
pressure  regions,  the  permeabilities  deviate  markedly  from  the 
otherwise  almost  linear  relation  between  log  D  and  l/T  as  ine 
temperature  decreases.  Figures  27  and  28  show  that  this  occurs, 
for  H2  and  D2,  only  in  filled  SNR  in  contrast  to  the  linear 


Fig.  25  -  Arrhenius  plot  of  permeability  isobars  for 
H2  and  unfilled  synthetic  natural  rubber. 
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Fig*  2?  -  Arrhenius  plot  of  permeability  isobars  for 
Ho  and  filled  synthetic  natural  rubber. 


Pig#  28  -  Arrhenius  plot  of  permeability  isobars 

for  Do  and  filled  synthetic  natural  rubber 


117 


plots  at  sufficiently  low  prossures.  These  deviations  are 
puzzling  and  possibly  experimental,  though  similar  observations 
were  made  by  Sobolev^O)  in  his  study  of  the  polyethylene- 
methyl  bromide  system, and  alike  trends  exist  in  our  S-P  curves, 
see  Figures  21  to  24, 

The  presence  of  crystallinity  causes  a  delay  of 
diffusion,  similar  to  that  of  caused  by  filler^0 »52)^  ^  number 
of  measurements  have  been  made  and  theories  have  been  postu¬ 
lated  based  on  the  general  physics  of  heterogeneous  systems 
containing  dispersed  particles^'^ ,  The  experiments 
showed  that  the  reduction  in  diffusion  was  a  function  of  the 
nature  of  the  gas,  so  that  a  simple  model  of  permeable  and 
impermeable  domains  is  not  sufficient (55 ) ,  The  tortuosity 
factor,  was  also  introduced,  but  still  the  agreement  with  the 
model  remained  unsatisfactory.  We  included  into  the  mechanism 
also  the  effects  of  the  presence  of  the  rigid  domains  on  the 
viscosity  of  the  metrix  and  thereby  obtained  an  improved 
agreement,  at  least  for  our  case.  However,  our  data  are  not 
good  enough  to  decide  whether  one  should  use  the  tortuosity 
factor  to  the  first  or  second  power,  as  also  debated  in  the 
literature ^6 ^  • 


B.  Activation  Energies  for  Diffusion  and  Permeation,  E,j  andET 


and  the  Heats  of  Solution,  AHs 


The  activation  energies  for  permeation  and  for 
diffusion  were  obtained  from  the  Arrhenius  equations# 


.118 


log  P  =  log  P0  - 
log  D  =  log  D0  - 


2.303  R  T 

Ed  1 
2.303  R  T 


at  various  pressures  for  the  four  samples.  They  are  listed 
in  Tables  3»  12,  21  and  30,  and  two  sample  plots  are  shown  in 
Figures  29  and  30.  Plots  of  the  activation  energies  against 
gas  pressures  are  shown  in  Figures  30  to  Jk,  It  is  seen  that 
the  activation  energies  are  independent  of  the  gas  pressure, 
although  the  values  of  Ep  for  filled  rubbers  scatter  and  may 
hide  minor  trends.  The  arithmetic  means  of  the  activation 
energies  for  various  gas  pressures,  Ep  and  E^,  for  all  samples 
are  listed  in  Table  1 6,  along  with  the  heats  of  solution,  *lHa, 
which  are  the  differences  of  these  mean  values. 

For  both  unfilled  rubbers,  the  activation  energy 
required  for  diffusion  of  D2  gas  is  about  350  cal/mole  greater 
than  that  of  for  H2  gas.  In  the  filled  rubbers,  the  case  is 
reversed,  and  the  activation  energy  for  diffusion  of*  K2  gas 
is  about  550  cal/mole  higher  than  D2  gas.  Although  these 
difference  are  small  on  an  absolute  scale,  the  least  square 
analysis  of  the  points  for  the  above  equations  (see  Appendix 
II)  demonstrates  the  reality  and  order  of  magnitude  of  the 
data. 

In  order  to  comment  on  these  findings,  it  should  be 
remembered  that  as  shown  in  Table  16,  the  heats  of  solution 
are  all  positive  except  for  hydrogen  in  the  filled  rubbers. 
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Pig*  31  -  Energy  of  activation  for  permeation  and  for 
diffusion,  Ep  and  E<«,  versus  gas  pressure 
for  H2  and  D2  in  unfilled  natural  rubber. 
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Pig*  32  -  Energy  of  activation  for  permeation  and  for 
diffusion,  Ep  and  E<j,  versus  gas  pressure 
for  H2  and  D2  in  unfilled  synthetic  natural 
rubber* 

A  denotes  D2  0  denotes  H2 


•NERGY  OF  ACTIVATION  (kcal  mole“ 


-r  — ■;■  yuc— ■w.tfcjijUBW 


fWr^jrafjmgjg!?«E|^tr-  '?£.' 3E 


123 


6  10  20  30  40  50  SO 

GAS  PRESSURE  ,  CM  OF  HG 


Fig*  33  -  Energy  of  activation  for  permeation  and  for 
diffusion,  Ep  and  E,»#  versus  gas  pressure 
for  H2  and  Djj>  in  filled  natural  rubber, 
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TABLE  16 


AVERAGE  ACTIVATION  ENERGIES*  FOR  PERMEATION  AND  DIFFUSION, 
Ep  AND  Ed,  AND  KEATS  OF  SOLUTION,  AHg 


Gas 

Unfilled 

NR 

Filled  NR 

AHS 

Ed 

EP 

Eq 

EP 

AHg 

5.36 

6,01 

0.65 

6.45 

6.73 

0.28 

H2 

5.03 

6.2? 

1.24 

7.03 

6.75 

-0.28 

6 

0.33 

-0.59 

-o.  58 

0.56 

Gas 

Unfilled 

SNR 

Filled  SNR 

Ed 

EP 

AHS 

Ed 

EP 

AHB 

d2 

5.32 

6.03 

0.71 

6,03 

6.38 

0.35 

h2 

4.93 

5.95 

1.02 

6.60 

6.26 

-0.34 

6 

0.39 

-0.31 

-0.57 

0.69 

*  Arithmetic  mean  in  units  of  kcal/mole. 

AHg  *=  Ep  ~  Ed 

6  is  the  difference  of  the  Ed’s  or  the  AHs's  for  D2  and  for 
H2»  i.e,,  the  difference  between  the  numbers  within  these 
columns. 

Standard  deviations  for  £d  and  Ep  are  listed  in  Appendix  II 
on  page  144. 


The  heat  of  solution  (mixing  or  absorption) ,  AHS,  in  the  un¬ 
filled  rubber  samples  includes  (see  the  discussion  above  on 
the  contributions  to  the  gas  solubility,  Srb»  in  rubber)  the 
heit  of  hole  formation  in  the  rubber,  AHy*,  the  heat  of  conden¬ 
sation  of  the  gases,  and  the  heat  of  internal  adsorption  of 
the  gas  to  the  rubber  chains,  £HC  andAHads  ,  i.e., 

AHS  «  AHh  +  AHC  +  AHads 

The  last  two  terms  in  the  equation  are  always  negative,  but 
the  middle  term  aHc  is  quite  small  in  the  case  of  these  gases. 
The  heat  of  mixing  of  gases  with  rubber,  AHa,  may  then  be  nega¬ 
tive  or  possitive  depending  on  the  heat  of  the  hole  formation* 
For  the  filled  rubbers,  there  are  two  additional 
terms  to  the  heat  of 'solution,  namely  ,  the  heat  of  adsorp¬ 
tion  on  the  filler  and  the  heat  of  gas  condensation  on,  or 
within,  the  filler  (absorption).  Therefore, 

AHS*  »AHS  +  AHads,  on  filler  +  ^cond,  inside  filler 

m  - 

The  last  two  terms  of  the  equation  are  strongly  negative. 

The  positive  AHS  listed  in  Table  16  implies  that  the 
energy  for  providing  the  volume  for  the  gas  in  the  rubber  is 
the  preponderent  term*  The  fact  that  the  heat  of  solution  for  Hg 
in  filled  rubbers  is  negative  comes  from  the  strong  adsorption 
of  H£  on  carbon  (£HadSl  on  fiHt;r9e”^Q0  cal/mole) .  On  the  other 

i 

hand,  the  small  positive  heat  of  solution  for  D-g  ^  filled 
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rubber  reflects  the  fact  that,  although  D2  is  also  more 
strongly  adsorbed  on  the  filler  than  on  the  rubber  chains,  the 
forces  are  not  quite  as  strong  as  for  H2t  in  other  words,  the 
potential  well  for  D2  on  the  carbon  is  less  deep.  As  a  result, 
deuterium  requires  a  lower  activation  energy  to  escape  from 
the  filler,  and  thus  from  the  filled  rubber  as  a  whole,  than 
H2  does.  The  energy  diagram  is  shown  in  Fig.  35* 

In  light  of  the  amounts  of  energy  involved,  about 
1500  cal/mole  both  penetrants  are  physically  adsorbed  on 
rubber  and  filler.  For  physical  adsorption  usually  only  a 
few  hundred  to  a  few  thousand  calories  per  mole  are  involved, 
due  to  the  operation  of  vender  Waals*  forces.  Chemi-sorption 
would  reveal  itself  by  much  stronger  binding  forces. 


Volume  of  Unfilled  Rubber 


Using  Frisch's  model,  as  described  in  Section  I  page 
22  «  the  average  dimension  of  the  unit  free  volume  or  free 
volume  increment  ,  V^,  in  the  rubber  can  be  estimated  by  the 
following  equation, 


18.8 

- - -  +  a 

(-6fHs  r 


Here,  6£Hg,  is  the  difference  of  the  heat  of  solution  for 
D2  and  H2  gas  in  the  same  medium,  i.e.,  6AHg  *  (AHg)^- 
and,  a  is  the  molecular  diameter  of  the  isotopic  species. 


The  quantum  mechanically  corrected  lennard-Joneo  potential 
diameter  of  hydrogen  (also  for  deuterium)  is  2.93  A^15^. 


t 


POTENTIAL 

ENERGY 


Fig.  35  -  Potential  energy  diagram  for  gas  on 
rubber  or  on  Carbon  Black 
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Using  the  average  value  of  AHS  listed  in  Table  16#  one  finds  ihat 
the  average  unit  free  volumes  for  the  unfilled  natural  and 
synthetic  rubber  are  about  51  and  64  f?  respectively  (see  . 
Table  1?). 

TABLE  1? 

AVERAGE  UNIT  FREE  VOLUME  OF  UNFILLED  NR  AND  SNR 


Sample 

( cal/mole ) 

.<ahs)H2 

(cal/mole) 

6AHS 

( cal/mole ) 

(X) 

Vf3 

(X3). 

NR 

650.0 

1,240 

-590.0 

3.7 

50.9 

SNR 

710.0 

1,020 

-310.0 

4.0 

63.9 

Thus,  the  average  unit  free  volume  in  NR  is  about  20#  smaller 
than  the  free  volume  in  the  SNR,  which  explains  the  higher 
permeation  rate  in  SNR  for  the  penetrants.  The  molecular 
basis  for  this  is  the  absence  of  crystallinity  of  SNR  which 
lacks  the  crystalline  nuclei  of  the  natural  rubber  latex  parti¬ 
cles,  i.e., SNR  is  truly  amorphous,  but  NR  can  be  considered 
as  lightly  filled  by  its  own  crystallites. 

Zf  one  attempts  to  calculate  the  free  volumes  of 
the  filled  rubbers  from  the  apparent  values  of  AHS,  this  is 
found  to  be  impossible.  The  reason  is  that  Frisch's  theory 
does  not  apply  to  this  case  in  that  the  theory  requires  6/Mis 
to  be  negative  in  order  to  get  a  positive  root  to  satisfy 
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the  equation  for  We  see  (Table  16)  that  for  the  filled 

rubber  $£HS  becomes  positive  which  in  a  literal  interpretation 
would  signify  that  H2  is  bonded  more  strongly  to  the  rubber 
than  D2f  a  physical  unreality.  This  dilemma  can  be  resolved 
by  assuming  that  the  adsorption  of  H2  onto  the  carbon,  where 
TT-complex  bonds  can  be  formed,  is  much  stronger  than  that  of 
D2.  If,  one  assumes  further  that  the  reversal  in  the  heats 
of  adsorption  is  entirely  due  to  the  presence  of  carbon  and 
that  the  solubilities  for  H2  end  02  remain  the  same  in  the 
matrix  of  the  filled  as  in  the  pure  rubber,  then  the  heats 
of  solution  in  rubber  and  thus  the  incremental  free  volumes 

1 

would  be  the  same  for  the  filled  and  unfilled  samples.  The 
intriguing  point  remains  that,  if  the  polarizability  difference 

between  H2  and  D2  played  a  part  in  the  adsorption  on  carbon, 

% 

and  if  this  were  true  for  the  heats  of  adsorption  on  the 
rubber,  then  the  heats  of  solutions  would  be  affected  not 
only  by  the  isotope  effect  but  also  by  binding  specificities j 
thus  the  outcome  of  the  free  volume  calculations  could  be 
changed.  Of  Course,  the  resulting  error  should  be  very  small. 

D.  Jump  Distance  of  Gas  in  Rubber 

The  jump  distance  of  the  gas  in  the  rubber  can  be 
calculated  by  Eyring*s  equation  (combined  from  equation  (33) 
and  the  Arrhenius  equation) , 


X  •  [t  2»1»AI  )*  D0  Vfl/3  ]  * 
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Using  the  average  values  of  D0  and  Vf  listed  in  Tables 
2,  5»  and  17#  the  culculated  jump  ,  distances,  X# 
for  the  two  gases  in  the  unfilled  rubber  are  shown  in  Table  18, 

TABLE  18 


JUMP  DISTANCES  FOR  HYDROGEN  AND  DEUTERIUM 
GAS  IN  UNFILLED  RUB3ERS  AT  40°C 


Rubber 

Gas 

D0  x  102 
(  cm  /sec  ) 

-1/3 

V 

(?) 

X 

(?) 

Hp 

5.49 

3.? 

21.1 

NR 

°2 

6.00 

3.7 

26.2 

h2 

4.95 

4.0 

•20.8 

SNR 

D2 

7.74 

4.0 

30.9 

The  jump  distances  should  be  larger  for  Dg  because  of  its 
greater  mass.  This  trend  is  further  reinforced  by  the  basic 
diffusion  constant,  D0,  which  is  also  larger  for  Dg.  As  a 
result,  the  jump  distances  of  D2  are  about  40  to  50#  higher 
than  those  of  Hg.  On  the  other  hand,  we  have  seen  in  Table 
14  that  the  difference  in  diffusivities  for  the  two  gases  in 
SNR  is  surprisingly  small.  This  cannot  be  explained  by  the 
larger  X  values  of  D2,  so  that  there  must  be  also  other  . 
factors.  For  example,  the  tortuosity  of  NR  due  to  its  crys¬ 
tallites  which  are  absent  in  SNR  may  affect  the  diffusivities 


of  gases  of  different  mass  differently.  Further,  if  D2  were 
markedly  less  adsorbed  within  the  non-crystallized  SNR  than 
it  is  in  NR,  while  the  adsorption  of  Hg  in  the  two  rubbers  is 
less  affected,  the  lack  of  difference  in  diffusivity  would 
become  more  understandable. 

E.  Conclusion 

Referring  to  the  Introduction,  this  investigation 
intended  to  find,  by  way  of  permeation  studies,  details  of 
structural  difference  between  NR  and  SNR  and  differences  re¬ 
sulting  from  filling  with  carbon  black.  In  particular,  the 
study  aimed  at  testing  and  applying  Frisch's  theory  and  at 
obtaining  further  information  on  -the  unit  size  of  the  free 
volume  in  rubbers. 

The  results  show  that  the  permeabilities  are  greater 
and  the  diffusion  is  faster  in  SNR  which  also  has  the  greater 
solubility  for  H2  and  D2.  This  agrees  with  the  notion^?) 
that  SNR  is  more  amorphous  than  NR. 

The  measurements  of  D  and  P  of  H2  and  D2  permitted 
then  to  calculate  the  free  volumes  and  jump  distances.  The 
results  are  physically  acceptable  and  in  line  with  the  data 
of  Frisch  and  Rogers^  1  \  and  Siegel  and  Eirich^0^.  Though 
the  differences  are  small,  SNR  shows  the  larger  values  for 
the  free  volume  and  jump  distances.  Since  the  trends  and 
consistency  of  the  data  support  their  statistical  significance 
(see  Ziegel) ,  the  findings  here  amount  again  to  a  confirmation 
of  Frisch's  theory. 
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The  results  obtained  on  the  filled  samples  permit 
interesting  conclusion  when  the  values  of  P,  D  and  S  are 
further  analyzed.  One  finds  that  one  can  relate  the  data  to 
those  of  the  unfilled  samples  if  one  takes  the  following  factors 
into  account t  the  blockage  of  part  of  the  cross  section  by  the 
filler,  the  imposed  tortuosity  of  the  diffusion  path,  and  the 
increase  in  viscosity  of  the  rubber  due  to  the  presence  of 
the  filler  particles.  Introducing  the  corresponding  correc¬ 
tions  brought  reasonable  agreement  between  experimental  results 
on  filled  and  unfilled  samples  even  though  the  corrections 
applied  can  not  be  tested  because  suitable  theories  have  not 
yet  been  developed. 

The  solubilities  for  H2  and  D2  were  found  to  be  low 
in  both  rubbers  and  the  same  is  true  for  the  adsorption  on 
the  carbon.  Thus,  the  later  is  probably  completely  wetted  by 
the  rubber.  As  a  consequence,  the  diffusion  of  H2  through 
filled  rubber  is  not  likely  to  be  modified  by  the  adsorption 
of  H2  on  the  carbon. 

Even  so,  gas  adsorption  on  the  carbon  is  about 
equal  to  the  absorption  in  the  rubber.  Therefore,  the  heats 
of  solution  are  more  negative  in  the  filled  rubbers.  Interes¬ 
tingly,  for  H2  is  much  more  negative  than  that  for  D2  so 
that  the  difference  in  overall  heat  of  solution  between  H2 
and  &2  “n  unfilled  rubber  is  not  only  wiped  out,  but  also 
becomes  positive  for  the  filled  rubbers.  Consequently,  the 
free  volume  increments  could  not  be  determined  in  the  filled 
samples,  but  it  is  reasonable  to  assume  that  the  free  volume 
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is  lowered  in  the  filled  rubber  since  the  diffusion  process  is 
impeded  beyond  reasonable  assumptions  about  a  diminished  cross 
section  and  increased  tortuosity* 

Comparison  of  the  heats  of  solution  between  the 
filled  and  unfilled  rubbers*  and  of  the  heats  of  adsorption 
on  the  carbon  derived  therefrom,  permits  to  estimate  the  heat 
of  adsorption  on  the  rubber*  Substracting  these  two  adsorp¬ 
tion  heats  from  the  overall  heat  of  solution  allows  to  estimate 
the  enthalpies  and  entropies  of  hole  formation  to  accomodate 
H2  gas.  This  energy  corresponds  to  that  of  separating  rubber 
molecules  during  a  swelling  process  where  H2  is  the  " solvent" • 

It  is  then  peculiar  that,  although  at  our  pressure  only  about 
one  K2  molecule  dissolves  per  several  thousands  cavities 
(which  explains  the  concentration  independent  of  the  diffusion) , 
the  heat  of  solution  is  positive  so  that  a  volume  expansion 
must  take  place  in  the  rubber*  This  could  be  explained 
by  assuming  that  the  existing  cavities  are  too  small  to 

accomodate  the  gas  molecules,  which  is  unlikely  in  view  of 

1/3 

the  Magnitude  of  Vf  ,  it  is  more  likely  that  the  intrusion 
of  K2  creates  some  changes  in  the  molecular  conformations 
which  are  akin  to  melting. 

Inasmuch  as  filled  rubber  has  a  heterogeneous  or 
composite  structure,  it  is  evident  that  the  method  applied 
here  could  also  be  used  more  generally  for  studies  of  compo¬ 
site  materials.  - 
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APPENDIX  I 

SAMPLE  CALCULATION  OF  PERMEABILITY,  P,  DIFFUSIVITY  D, 

AND  SOLUBILITY,  S 

The  application  of  equations  (42)  and  (43)  to  calculate 
P  and  D  is  illustrated  in  the  following  example.  Let  us  examine 
the  versus  t  plot  as  shown  in  Pig,  12  for  the  permeation  of 
D2  gas  at  70°C  through  unfilled  natural  rubber  of  thickness 
0,2101  cm  under  a  gas  pressure  of  17,60  cm  of  Hg,  The  slope 
of  the  linear  portion  of  the  plot  is 

b  *  3,2786  micron/min, 

and  the  value  of  the  intercept,  MaM,  is  calculated  from  the 

slope,  b,  and  a  value  of  tk  with  its  corresponding  value  of 

« 

pressure  pk  as  indicated  in  the  figure.  Therefore, 

a  *  pk  -  btk  «?  60  -  3,2786  x  23*4 
*  60  -  76.72 
*  -I6.72 

Substituting  this  value  together  with 

P  =  0.02  micron  (initially  measured  residual 

pressure  ) 

L  =  0.2101  cm 


into  equation  (42),  one  obtains, 


^*.-500.  \-fm  ~ 


rsizx  irtr-y^isl  3E£&E&r=s:  i?ax^T 
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L2  x  b  (0,21 01. ) 2  x  3.278* 

*  360  (P0-  a  )  *  3*o[o.o*  -  (-l*.?2)] 

*  2.4o  x  10"-  nm2/p<?c» 

Th*  slope  b  in  *5g.  12  als"  represents  the  steady 
state  flux  Qs  in  units  of  micron  per  minute.  Assuming  ideal 
gas  behaviour*  the  number  of  moles  of  gas  transported  per 
second  is 

b  x  10“4  Vi 

q  .«  — - -  x  —  moles/sec 

8  76  x  60  RT 


where  T  is  the  temperature  of  the  gas  collector.  Q_  is  custo- 
marily  expressed  in  units  of  standard  cc  of  gas  defined  as 
cc  of  gas  at  25°C  and  one  atmosphere  pressure.  Therefore, 


=  b  x  10"^  .  Vj  x  2?B 

8  76  x  60  X  T 


b  x  10~^x298  x  Vt 
76  x  60  x  T 


Std.  cc/sec 


Substituting  this  value  of  Qg  into  equation  (43),  one  obtains 
the  permeability, 

L  b  x  lO^x  Vi  290 

P  *  (  Ap/?6  )  x  A  X  76  x  60  X  T 

*  8.51  x  10"?  Std.  cc  gas  sec"1 cm"1 atm"1 


1.1*0 

where  Ap  s  1?.6  cm  Hg 

L  =  0.2101  cm 
b  =  3*2?86  micron/min 
Vx  =  37*9  cc 
A  =  2.85  cm2 

T  =  29S°K  (when  the  collector  is  at  the  room 
temperature  ) 

Finally,  according  to  equation  (11),  the  solubility, 

S  =  P/D  -  3./S5  eft  gas  em"3  atm-1 


APPENDIX  II 


l4l 


A.  SAMPLE  CALCULATION  OF  ENERGIES  OF  ACTIVATION  FOR- 
PERMEATION  AND  DIFFUSION(Ep  and  Ed) ,  AND  HEAT  OF  SOLUTION(AHg) 

BY  THE  LEAST  SQUARE  METHOD 

B.  STANDARD  DEVIATION  FOR  E?.  AND  E 

d  P 

A.  The  logarithmic  form  of  Arrhenius  equation  for  diffu¬ 

sion  and  for  permeation  is  as  the  follows  (quoted  from  P.  53) 


and 


log  D  *  log  D0 
log  P  *  log  P0 


Ed  1 

2*303  R  T 

Ed  1 

2*303  R  T 


Both  above  equations  can  be  simply  represented  by  a  straight  line 
equation  as  * 

y  *  a  +  bx  (46) 

where  y  *  log  D  or  log  P 

a  *  log  D0  or  log  P0 
b  *  Ed/2.303  R  or  Ep/2*303  R 
x  *  i/r 

♦ 

For  a  set  of  experimental  data,  the  deviation  of  each 
datum  point,  from  a  best  straight  line,  <^,  can  be  expressed  by 


«i  •  n  -  a  -  bx* 


Summation  of  the  square 


2 

( 


of  each  deviation, 
n 


i=l 


obtains 

2 

bxi) 


* 
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i 

l- 


l 

S 


K 

i 

K 

| 

£ 


E 


| 


li 


RT. 


2  2 

=  (y1  -  a  -  bx1)  +  (y2  -  a  -  bx2)  +  . 

*  (yx2  +  y22+“*)  +  b2^2  +  x22  +  ”••) 

-  2a  {y1  +  y2  +  •»•)  -  2b  (x1y1  +  x2y2  +  •  ••) 
+  2ab  (x  +  x2  +  •••)  +  na 

minimize  <2  with  respect  to  "a"  and  Mb"  respectively,  giving 
ZL  m  -2  (y1+y1  +  ...)  +  2b  (  xx+x2  +  ...)  +  2na  *  0 


or 


na  «  b  e  xi  «  Eyi 


(47) 


2 

~~  =  2b  (  x12  +  x22  +  ..*)  -  2  (x^  +x2y2  +  ...) 


+  2a  (  +  x2  +  •  • ,)  «  0 


or 


n 


n 


n 


Ixi+blx 


i2  *  XXiyi 


(48) 


i=l  i=l  i=l 

Solving  for  Ma"  and  "bM  from  above  two  equations,  results 


n 

y  2-*i  ■  *  ]£  Xiyi 

iel  i*l 


n  n 

Zxi2-xZ  xi 


(49) 


i-1 


i=l 
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and  hence,  the  slope,  "b",  and  the  intercept,  Ma"  of  the 
straight  line  can  be  separately  calculated  by  equations  (49) 
and  (50)  for  a  set  of  data.  Further  according  to  equation  (46), 
values  of  DQ,  Ed  and  Ep  i.re  readily  calculated  for  each  gas 
pressure. 

To  illustrate  this,  the  permeation  data  of  Table  2 
on  P.  55  for  H2  gas  transport  through  the  unfilled  natural 
rubber  at  gas  pressure  of  1,1  cm  Hg  are  quoted  to  calculate 
E(j,  D0  and  Ep  as  the  following* 


Temp.- 

(  °0) 

D  x  105 

(cm2/sec) 

< 

H 

11 

•H 

X 

Yi  *  log  D 

40 

1.49 

3.2  X  10-3 

-4.8269 

55 

2.54 

3.05x10*3 

-4.5960 

70 

3.18 

2,92  x 10"3 

-4.4981 

85 

4,43 

2.79 x  10-3 

-4.3537 
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Applying  equations  (49)#  (50)  and  (46),  one  obtains, 
a  **  log  D0  *  -1.148 

or  D0*  7»1  x  10“2  cm2/sec- 

b  «  -1.145  x  103 

and  Ed  *  5*24  kcal/mole 

Similar  calculation  for  energy  of  activation  for  permeation 
Ep  arrives  at 

E  *  6.34  kcal/mole 
P 

Finally,  the  heat  of  solution,  AH8  is 

AHS  "  Ep  -  Ed  ®  1 1 10  kcflil/nioX6 

B.  STANDARD  DEVIATIONS  FOR  Ed  AND  Ep 

_ Unfilled  NR 

Ed  Ep 

D  2  5.36*0.21  6.01  to. 19 
H2  5.03*0.17  6,27  i  0.06 


_ Filled  NR 

Ed  Ep 

6,45+0.17  6,73  to, 20 

7.03*0,15  6,75*0.28 
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Unfilled  SNR _  _ Pilled  SNR 


*d 

Si 

h 

d2 

5.32  10.22 

6.03  *0.20 

6.03  ± 0.19 

6.38*0.30 

H2 

4.93*0.10 

5.95*0.13 

6.60  t o.l6 

6.26  1 0.26 

amuiiHW-r,! 
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